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Abstract 

Recently  Kralchnan  (1959)  has  propounded  a  theory  of  homogeneous 
turbulence,  based  on  a  novel  perturbation  method,  that  leads  to  closed 
equations  for  the  velocity  covariance.   In  the  present  paper,  this  method 
is  applied  to  the  theory  of  turbulent  diffusion  and  closed  equations  are 
derived  for  the  probability  distributions  of  the  positions  of  marked 
fluid  elements  released  in  a  turbiilent  flow. 

Two  topics  are  discussed  in  detail.   The  first  is  the  probability 

distribution,  at  time  t,  of  the  displacement  of  an  element  from  its 

initial  position.   In  homogeneous  flows,  this  distribution  is  found  to 

resemble  that  for  classical  diffusion  but  with  a  variable  coefficient  of 

diffusion  which  is  proportional  to  v  t  for  t  «  t/v  and  which  approaches 

a  constant  value  =  Iv     for  t  »  l/v        (^  =  macroscale,  v  =  r.m.s.  turbu- 

o  '  o  o 

lent  velocity) . 

The  second  topic  treated  is  the  joint  probability  distribution  of 
the  displacements  of  two  fluid  elements.   Particular  attention  is  focused 
upon  the  probability  distribution  of  relative  displacement,  i.e.  lUchard- 
son's  distance -neighbour  function.   This  is  found  to  be  Gaussian  for 
separations  r  which  are  large  (»  t) .      For  smaller  separations  (r  «  l) , 
its  behavior  at  high  Reynolds  numbers  is  foxmd  to  be  quite  well  expressed 
in  terms  of  a  variable  diffusion  coefficient  K(r,t),  as  suggested  by 
Richardson  (1926).   For  all  but  extremely  short  times,  K(r,t)  is  found 

to  depend  only  on  r  and  on  the  form  of  the  Inertial  range  spectrum  E(k). 

2     -5/2 
On  assuming  E(k) oC  v  t(kt)     as  results  from  Kraichnan's  approxima- 

tion  (1959),  one  finds  Kir)  oC  v  l{r/l)-^'    .   On  the  basis  of  similarity 


ii 


2  -'5/3 

arguments   of  the  Kolmogorov  type^    -which  give  E(k)  oC  v     l{'kt)    ^'    ,    one 

i|./5 
finds  K(r)  oC  v  1{t/1)  as,    in  fact,    Richardson  originally  proposed. 

The  dispersion    \t     /is  proportional  to  I   (v  t/t)      on  Kraichnan's  theory^ 

while    \r    /   oC   t   {v  t/l)      on  the   similarity  theory.      This  illustrates 

that  the  behavior  of    <^r    ^  is  very  sensitive  to  the   spectrum. 
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1.    INTRODUCTION 

The  aim  of  the  theory  of  turbulent  diffusion  is  to  determine  in  a 
statistical  sense  the  migration  of  marked  particles  as  they  are  carried 
along  with  a  turbulent  flow.  Like  molecular  motion  in  a  dilute  gas  of 
discrete  particles,  turbulent  diffusion  is  a  linear  process  if  the  con- 
vected  particles  have  no  reaction  on  the  flow;  i.e.  the  probability 
distribution  for  the  position  of  a  marked  particle  in  space  obeys  the 
superposition  rule  and  changes  in  time  according  to  a  linear  equation. 
Unlike  classical  molecular  motion,  the  motion  of  neighboring  fluid  ele- 
ments in  a  continuiim  is  correlated,  although  one  expects  that  over  dis- 
tances large  compared  to  the  macros cale  I   cf  the  turbulence  this  corre- 
lation is  weak  and  that  elements  separated  by  such  distances  move  almost 
independently.   Furthermore,  unlike  classical  molecular  diffusion, 
turbulent  diffusion  is  not  a  Markoff  process.  However,  one  expects  that 
over  times  large  compared  to  l/v   ,  where  v  is  the  root  mean  square  fliild 
velocity,  the  fliiid  elements  will  suffer  many  essentially  uncorrelated 
deflections  by  the  energy  containing  eddies  and  that  accordingly  their 
motion  over  such  long  times  will  be  almost  Brownlan.   Under  such  circum- 
stances, one  expects  that  the  spread  of  marked  particles  carried  by  the 
fluid  will  indeed  resemble  classical  diffusion,  and  that  it  will  be 
possible  to  define  a  coefficient  of  eddy  diffusivlty.   An  analytical 
basis  for  these  qualitative  observations  is  given  in  Sec.  2. 

Particles  which  start  out  simultaneously  at  nearby  points  have 
closely  similar  histories  in  any  one  realization  of  the  turbulent  flow 
(and  over  times  which  are  not  too  long).   For  such  times,  their  relative 
motion  is  unaffected  by  eddies  whose  spatial  scale  is  large  compared  to 


the  Initial  separation;  such  eddies  give  nearly  equal  displacements  to 
the  tvo  points.   The  change  in  the  separation  r  is  governed  by  the  smaller 
eddies,  particularly  those  vhose  length  scale  is  of  the  same  order  as  r. 
Thus,  in  a  flow  of  high  Reynolds  number,  we  expect  that  while  the  particles 
are  separated  by  a  distance  appropriately  small  compared  to  I   their  rela- 
tive diffusion  will  be  governed  by  the  inertial  range  spectrum  of  the  tur- 
bxilent  flow,  and  will  be  unaffected  by  the  structure  of  the  energy  contain- 
ing eddies.   An  analytical  basis  for  these  surmises  is  given  in  Sec.  5, 
and  a  form  is  proposed  for  the  variable  diffusion  coefficient  introduced 
by  Richardson  (1926). 

2.    DIFFUSION  FROM  A  FIXED  SOURCE 

2.1  Methods  of  approach 

There  are  two  main  ways  of  attempting  to  give  an  analytical  frame- 
work to  the  qualitative  arguments  of  Sec.  1.   In  the  Lagrangiaa  framework 
(as  distinct  from  the  Eulerian  approach  to  be  described  presently),  pro- 
bability distributions  are  defined  for  the  displacements,  velocities,  etc., 
of  given  marked  particles,  and  the  relationships  between  them  are  studied. 
For  exemple,  let  G(x,t|x  ,t  )dx  be  the  probability  that  a  fluid  particle 
lying  at  the  point  x  at  time  t  should,  at  the  later  time  t,  lie  within 
a  voltmie  dx  at  the  point  x.   Let  V(x,u,t|x  ,t  )dxdu  be  the  probability 
that  this  same  particle  should  at  that  time  lie  within  dx  and  have  a 
velocity  between  u  and  u  +  du.   Then  it  is  not  difficult  to  show  that 

J-  G(x,tlx  ,t  )  +  ^   /  V(x,u,t|x  ,t  )u.du  =  0.  (2.1) 

dt  ^-.'  '~o  o    dx  J        ^~.'-^'    '-0  o  1  ~ 

(in  2.1  and  elsewhere,    we  use  the   sintmiation  convention.) 


5. 


Through  the  hydrodynamical  equations,  it  is  possible  to  derive  a  similar, 
though  far  more  involved  equation,  relating  V  to  another  probability 
function,  and  this,  in  its  turn,  to  yet  another.   It  woiild  be  necessary 
to  close  this  hierarchy  of  equations  in  some  vay  in  order  to  obtain  from 
it  an  evaluation  of  G(x,t|x  ,t  ). 

The  second  approach  (advocated  by  G.K.  Batchelor,  1952b;  see  also 
W.H.  Raid,  1955;  P-H.  Roberts,  1957)  Involves  reformulating  the  problem 
of  finding  Lagrangian  probability  functions  as  one  of  determining  certain 
Eulerian  moments.   A  passive  scalar  quantity  \|/(x, t)  is  introduced  vhich 
satisfies  the  equation 

1^  t(x,t)  +  ^  0|'(x,t)u.(x,t)]  =  0,  (2.2) 

1    ~ 

where  u. (x,t)  is  the  velocity  field.   This  quantity  is  therefore  carried 

1  r— 

by  the  turbulent  fluid,  but  does  not  affect  its  motion.   It  is  clear  that 
if  we  take 


t(x,t^)  =  S(x  -  x^),  (2.3) 

where  6(x)    Is   the   three   dimensional  Dirac  5-function,    then 

)|/(x,t)    =  6(x  -   x^),  (2.1|) 

where  x  is  the  position  at  time  t  (in  this  particular  realization  of  u) 
of  the  fluid  element  which  was  initially  at  x  .   By  averaging  over  all 
realizations,  we  see  that  the  solution  of 

^<^(x,t)>  +  J-<^(x,t)u.(x,t)>  =0,  (2.5) 


vMch  satisfies  the  initial   condition 

<^(x,t^)>    =  6(x  -  xj,  (2,6) 

Is 

<;i|r(x,t)>    =  G(x,t|x^,t^).  (2.7) 

Through  a  similar  though  more  complicated,  equation,  <^\|m.>  is  related 

to  higher  moments,  and  so  on.  Again,  we  have  an  hierarchy  of  equations 

vhich  must  be  closed  in  some  way  in  order  to  evaluate  G(x, t|x  ,t  ). 

In  the  present  paper,  we  adopt  the  second  of  these  two  approaches 

and  employ  an  approximation  method  for  closing  the  equations  devised  by 

Kraichnan  (1959).   However,  before  presenting  this  approximate  analysis, 

we  shall  derive  some  results  which  eire  asymptotically  exact  for  small 

t-t  . 
o 

2.2  Exact  results  for  short  times 

For  t-t  «  t/v  ,  the  fluid  particles  are  simply  swept  from  thfeir 
o       o 

points  of  origin  with  whatever  velocity  the  -turbiilent  fluid  happens  to 
have  at  the  moment  of  their  release;  i.e. 

V(x,u,t|x^,t^)  =  P[u(x^,t^)]6[(x-x^)  -  u(t-tj],  (2.8) 

where  P[u(x  ,t   )]    is  the  probability  distribution  function   (p.d.f.)    of 
u  at  position  x     and  time  t   .      Thus,   by  (2.1),    or  by  inspection 

m,  ~.0  O 

G(x,tlx  ,t    )    =  — i 5P(|-|^).  (2.9) 

**        -^     °  (t-tj5       \t-t^j 


This   res\ilt  can  also  be  deduced  from  the   formal  solution 


t( 


.|df   ^    ]u^(x,f)|    at"  ^   [u.(x,t")^(x,t^)] 


(2.10) 


vhlch  one  obtains  from  (2.2)  by  integration  and  iteration.  When 
t-t  «  t /v  ,  then  u(x, t)  =  u(x  ,t  ).  and  it  follows  that 


t(x,t)= 


l-(t-t  )u^(x  ,t   )^  +  -^(t-t  )  u,  (x  ,t  u.(x  ,t    )  ^ — ^ — 

^   o^  i  ~o'  o'^dx.    2i^   o'   i^-o'  o  j^~o'  o'  dx.dx. 

—  1  '^         1   J 


ij/(x,t  ) 


(2.11) 


Hence,  by  (2.6)  and  (2.7), 


G(x,t|x  ,t    )= 

^«'   'wO   o 


l-(t-t  )  <u,  (x  ,t  )>  ^ 
^    o   \  l^-O   O  ^  OJ 


IP  ■\'- 

+  ■?rr(t-t  )  'vu.(x  ,t   )u.(x  ,t  )/  ^ — ^ .... 

2:^    o   \  i^~o  o  J  "O  o    dx_,ax. 

i   J      - 


6(x-x  ) 


(2.12) 


That  this  is  equivalent  to  (2.9)  can  most  easily  be  seen  by  expressing 
the  result  in  vave-vector  space,  writing 


G(x,tlx  ,t  )  =  /  G(k,t|x  ,t    )e 

^~'      '-O   O     J     ^-'   '-0   o 


ik.(x-x  ) 


dk. 


(2.13) 


6. 


Then  (2.12)  is  equivalent  to 


G(k,t|x  ,t    )  =  — ~   l-ik.(t-t  )  /u,  (x  ,t    )> 
^^'    i»o  o    /p  kI—   ■""    °     i  "O  o 

1  2 

-  -jr-p  k.k.(t-t  )  /u.(x  ,t    )u.(x  ,t   )  y  +   ... 
2:   i  j^    o   \  1^-0  o  .I  ~"o  o  ' 


(2.11^) 


(2.15) 


where 


P(T,)  =  rP[u(x^,tj]e  '-'-du. 


(2ol6) 


is  the  characteristic  fimction  for  the  distribution  of  velocity  at  x 

■'    ^o 

and  t  .  Equation  (2.9)  is  simply  the  inverse  of  (2.15). 

This  second  method  of  establishing  the  behavior  of  G  at  small  times 
brings  out  some  noteworthy  features.   If  (2.12)  is  cut  off  after  any 
finite  n-umber  of  terms ^  it  implies  that  G(x,t|x  ,t    )  vanishes  identically 
for  non-zero  x-x  .   On  the  other  hand,  if  (2.1^^-)  is  cut  off  after  a 

—  -O  7      \      / 

finite  number  of  terms,  the  resulting  expression  for  G  diverges  for 
large  k,  or  large  t-t  .   One  concludes  that  any  reasonable  approximate 

kiv  O 

solution  for  the  full  space-function  G(x, t|x  ,t  )  must  include  terms  of 
orders  of  the  formal  expansion  (2.10).   Even  for  very  short  times  the 
formal  expansion  is  only  useful  because  we  happen  to  be  able  to  sum  it 
to  all  orders . 


7. 


Hovever,  if  the  moments 


<^i>  =/ 


Ax^Gfx.tlx  ,t    )dx  =  i(2rt)- 


^  G(k,t|x  ,t  ) 
ok.   ^~  '-0  o  , 

>—    1  -•  K=0 


< 


Ax,AXj> 


Ax,Ax,G(x,t|x  ,t    )dx  =  -(2rt)' 

1    J   —    '  —  O   O   «" 


"^k^^(^^Lv^o) 


i-  i   J 


-ik=o 


(Ax  =  X  -  X  )  are  expanded  by  means  of  the  formal  solution^  the  resulting 

•»         «MK         <M«0 

series  appear  to  converge  for  all  t-t  ,    although  the  convergence  is  poor 

unless  t-t  «  t/v   .      Equation  (2.1^)  shows  that  for  small  t-t  , 
o     '  o  ^     '  o 


/Z^x.Ax.N  =  U.  .(x  ;t  ;x  ,t  )(t-t  )' 
^   i   j/     ij  u,o  o'..o  o'^    o 


(2.17) 


where 


U.  .(x,t;x  ,t    )   =  <'u^  (x,t)u.(x  ,t  )'> 
ij^fc-'  '.^o'    o  \  iotas'  '    .1  -o'  o  '^ 


is  the  velocity  covariance.  In  the  Isotropic  case^  therefore. 


<'Z\x.Ax.\  =  v,  5,  .(t-t  )^, 


(2.18) 


where  v   is  the  mean  square  of  any  component  of  velocity  at  position  x 


and  time  t  . 
o 


There  is  fairly  strong  experimental  evidence  (see,  for  example, 
Batchelor,  1953?  Ch.  8)  that  P  almost  always  is  closely  Gaujssian.   In 
this  case,  G(x, tlx  ,t  )  must  be  closely  Gaussian  also  for  short  times. 

«•    'mO   O 


8. 


Then  by   (2.1^1)    and   (2.9)   we  have 


G(k,tlx    ,t    )    =  — ^  exp     -  i  U     (x  ,t    ;x    ,t    )k  k   (t-t    r 


(2«) 


and 


(2.19) 


Q(^^li^o-^o^"..    ,3/2.      ^    ,5, 


(2jt)^^    (t-t^)-^(det  U^    ) 


V2 


X     exp 


1  -2 
7c  u.  ,  A  X.    A  x.(t-t    ) 

2  ij  1  J  o      _ 


(2.20) 


where  u.  ,   Is   the   cof actor  of  U.  .(x    ,t    ix    ,t    )    and  det  U_,  .    denotes   the 
ij  ij^-o     o  -o     o  Ij 

determinant   of  these   quantities.      In  the   isotropic    case. 


G(k,t|x   ,t    ) 

n>  bkO       O 


(2rt)- 


exp 


1  1  2     2,,    ^    ,2 
-  -  k  v^   (t-t^) 


and 


(2.21) 


G(x,tlx   ,t    ) 


[2:tv/(t-tj2] 


exp 


-(x-x   )^-/2v,^(t-t    f 

^m,     VO  1       ^  O 


(2.22) 


2.5     An  integro-differentlal  equation  for  G(x, tlx   ,t    ) 

' %r  VfO  O 

When  t-t   is  not  small  compared  to  t/v  ,  the  approximation  (2.11)  to 
the  formal  solution  (2.10)  is  Invalid  since  it  is  no  longer  legitimate  to 
Ignore  the  space  and  time  variation  of  u.(x,t).   It  is  nevertheless  possible 
to  effect  a  partial  summation  of  (2.1)  which  includes  terms  from  every  order 
in  the  expansion^  and  is  such  that  the  resulting  expression  for  G  converges 


for  large  k.      The  integral  equation  for  this  approximate   form  of  G  is 
derived  in  Appendix  A  in  tvo  ways.      The   first  makes  use   of  Kralchnan's 
direct  interaction   approximation   (Kraichnan,    1959)*      In  the   second  method, 
the   same  result  is   derived  by  discarding  or  retaining  terras   in  the   formal 
expansion   according  to   a  certain  selection  criterion.      The  terms   retained 
are  of  all  order.      Each  of  these  methods   supposes   that  the  velocity  field 
is   spatially  homogeneo\is,   but  this   restriction   can  be   removed  by  the 
application  of  a  more  general  method  due  to  Kraichnan  (1960).      The   final 
result,    for  the   case  vhere  the  mean  field   <'u(x,t))>     vanishes,    is  the 
equation 


|^G(x,tlx^,tJ   =    f    dt'y^dx'U..(x,t;x^,t^) 


>^     ^  G(x,tlx',f)  ^  G(xStlx^,tJ.  (2.25) 

i  J 


Here,  as  later,  the  fluid  velocity  is  supposed  incompressible: 

Su.(x,t)/ax.  =  0.  (2.21^) 

X  «»     1 

(The  compressible  case  can  be  treated  by  similar  methods.)   In  the  case 
of  statistically  homogeneous  and  stationary  flovs,  we  may  write 


^(^^l^o'^o^  =  ^^l-^o'^-\^' 


U.  .(x,t:x  ,t  )  =  U.  .(x-x  ,t-t  ), 


(2.25) 


10. 


and^  upon  a  partial  integration^  (2.23)  becomes 


^=        * 


G(x,t)  =  V  ^s     /   dt'  /  dx'U.  .(xSt')G(xSt')G(x-xSt-t').      (2.26) 
■^   o 


Multiplying   (2.26)  by  x.x     and  Integrating  the   right-hand  side  by 
parts  ^   we   find 


^<x.x.>=    /x,x.|G(x.t)d.=   2    f    df|dx'U..(x',f)G(x',f); 

o 

i.e. 

<x.x.>    =   2tK^^(t),  (2.27) 

where  the  variable   "eddy  diff\isivity"   tensor  K.  .(t)   is   defined  by 

J-  O 

t 
K.  .(t)  =  X  f    dt'(t-t')  /  dx'U.  .(xSt')G(xSt').  (2.28) 


For  short  times  (t  «  l/v   ),    (2=27)  and  (2.28)  agree  with  (2.17),  if 
we  assume  G(x',t')  is  negligible  unless  |x'|  «  I.      Then, 


K.  .  =  i  t  U.  .(0,0),     t  — >  0.  (2.29) 


For  large  times  (t  »  l/v   ), 


11. 


. ,=  /  dt'  /  da'U. .(x',t')G(x',t' )  =   ^  ^  /  dt'  /  dkU. .(k,t')G(-k,t'). 


'   '    (2«) 

o  ^ 


(2.50) 


Let  us  assume  that  the  diffusion  for  time  >  t/v  is  dominated  by  the 

'^     o 

energy- containing  eddies.   Then  it  is  reasonable  to  suppose  that  the  integral 
in  (2.50)  should  depend  only  on  the  parameters  I   and  v  ,   whence,  by  dimen- 
sional reasoning,  we  must  have 


K.  .   — >  constant  of  order  tv  ,   t  — >  °° 
ij  o' 


In  the  Isotropic  case,  we  have  (setting  x  =  |x| 


IJ      IJ  5 


j      dt'(t-t')  /  dx'x'V^(x',t')G(x',t').       (2.51) 


Equations  (2.27)  and  (2.28),  and  their  generalizations  for  flow  which 
are  not  steady  or  homogeneous,  imply  that  the  Lagrangian  correlation  func- 
tion <u.  (t)u.(t  )N    (cf„  Taylor,  1921)  is 
1     J   o  ''  L 

<u.(t)u.(t  )>^  =   /  dxU.  .(x,t;x  ,t  )G(x,t|x  ,t  ).  (2.52) 

^  1^    j^  o  ^L   J      -   ij  ~-   —o  o  V  'wo  o  \      ^   I 

This  should  be  compared  to  the  exact  result 

<u.(t)u.(t  )>_  =   /  dx<u.(x,t)u.(x  ,t  )ii/(x,t|x  ,t  )>  .  (2.55) 
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where  \(r(x,t|x  ,t    )  denotes  the  unaveraged  Green's  function;,  i.e.  the 
N*   '  **o   o 

solution  of  (2.2)  satisfying  (2.3).   Equation  (2.52)  clearly  may  be 
obtained  from  (2.53)  on  the  assumption  that  the  \jf  and  u  fields  are  uncor- 
related  (see  also  Corrsln,  I960)  .   However,  if  one  decides  to  use  an 
assumption  of  this  kind,  the  results  depend  very  much  on  what  stage  of 
the  analysis  is  chosen  for  its  application.   For  example,  had  we  assumed 
that  ijf  and  u  are  uncorrelated  in  (2.2),  we  would  have  obtained  the  absurd 
result  d  <\|/(x,t)^  /St  =  0.  Moreover,  if  the  approximation 

<u,(x,t)u.(x^,t  )ij/(x,t|x^,t  )>  =  U  (x,t;x  ,t  )G(x,t|x  ,t  )         (2.3^) 

1  •-       J  "..O    O     •»»     ukO    O  Ij  m-    "uO    O     M^  "«.0    u 

is  used  to  close  the  hierarchy  of  equations  which  arise  from  (2.2)  and 
(2.5),  the  resiat  is  not  consistent  with  (2.52)  (cf.  Roberts,  1957)-   Thus^ 
the  direct-interaction  approximation  is  consistent  with  (2.5^)  only  if 
(2.5^)  is  used  in  a  particular  way. 

2.i<-   Solution  of  the  integro-differential  equation  for  short  times 

For  t  «  -t/v  ,  we  may  assume  that  the  G  factors  in  the  integrand  of 
(2.25)  are  negligible  ixnless  x-x  and  x'-x   are  both  small  compared  with  I. 

^  ^'  00  iwi    -O  <*v      wO 

Thus,    n    .(x',t',x   ,t    )    in  the   integrand  may  be   replaced  by     U      (x   ,t   ,x   ,t    ) 
'      ij^v.         '«^o     o  o  ^  Ij  ».o     o  kfcO     o 

2 
or,    in  the   isotropic  case,   by  v     5. ..      It  follows   that 


t 
^  G(x,t)   =  v^^  V^      r    df    rdj^'G(x',t')G(x-x',t-t').  (2.35) 

o 


Saffman   (\inpublished)    has   used  the   approximate   result    (2.32)    to   derive   an 
estimate   for  K   (see   Sec.    2.5  below). 
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Taking  the  Fourier  transform 


G(x,t)  =  /   G(k,t')e^^'-dk  , 


we  find 


^"-^t 


G(k,t)  =  -(2n)\^\^   /  dt'G(k,t')G(k,t-t') 


(2.56) 


By  taking  the  Laplace  transform  of  (2.56),  or  by  using  a  result  of  Watson 
(19^^,  Sec.  12.2,  eq.  5),  we  find 


1    Ji(2kv^)         ,  , 
G(k,t)  =  — ^   \     /   ,   (k  =  k  ) 
(2n)5    l^V 


(2.57) 


Inverting  by  using  a  further  result  of  Watson  (19'<-^^  Sec.  13.^2,  eq,  (^)),  we 
find 


/> 


(atv^t)' 


1     2^2      2I 

hr     t  -  X 


-1/2 


,  if  X  <  2v  t, 


G(x,t)  =  < 


(2.58) 


V 


,  if  X  >  2v  t. 


The   corresponding  probability  distribution  D(x  ,t)    for  displacement   along 
the  X-axis    (whose   direction  may  be    chosen   arbitrarily)    is   related  to 
G(x,t)   by 


SD(x,t)/ax  =   -2rtxG(x,t), 


(2.59) 


Ik. 


and  so,  in  the  present  case. 


/ 


2it(v^t)   L 


I   2^2     2 
W^t   -x^ 


1/2  , 


if  x^  <  2v^t, 


D(x^,t)  =   < 


,if  x^  >  2v  t 


(2.1^0) 


Both  G(x,t)    and  D(x   ,t)    are   everywhere  non-negative,    but  they  are   zero  beyond 

a  distance   of  2v  t   from  the   source.      That   this  behavior  is   not   restricted  to 

3/2  - 
small  times    can  be   seen  in   a  rough  way  from  the   fact  that  k  G(k, t)    oscil- 

lates with  wavelength  2v  t   as  k  — >  «>    .      A  more    convincing  demonstration   can 
be   given  directly  by  induction  from   (2.26).      One   supposes   that  G(x',t') 
vanishes  when    \x' \    >  2v  t'    for  all  times   t'    ^   t"    and,    by  approximating  to 
(2.26)    for  points   near  the  wave-front    Ixl    =   2v^t"    ,    one    can  show  that  this 
implies   that   G(x,t)    vanishes   for    Ixl    >  2v  t,    where   t   is   infinitesimally  great- 
er  than  t" .      The   finite  maximum  "propagation  speed"    exhibited  by   (2.38)    there- 
fore persists   for  all  times.      This   is    also   consistent  with   (2.9)   which  implies, 
in   the  present    case,    that  the  p.d'. f.    of  velocity  has    a  sharp   cut-off  at 
u    =   2v    .      These  properties    are  physically  unreasonable.      They  arise   in  some 
way  from  the   direct   interaction   approximation. 

A  related  unrealistic   feature   of   (2.56)    is   that  G(x,t)   — >  <»   as 
X  — >  2v  t   -   0.      However,    the   singularity  is  weak,    since 


2v  t  +  0 


/ 


4rtx^G(x,t)dx  =     -  (2  -  -^) 
\     '     '  jt     \  V    t 


1/2 


->  2v^t    -    0) 
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i.e.  the  singular  outgoing  wave  front  does  not  even  carry  a  finite  inte- 
grated probability.   It  is  perhaps  worth  mentioning  that  this  singular 
type  of  behavior  is  not  unknown  in  random  walk:  problems  although  it  is 
generally  smoothed  away  after  a  few  steps.   (For  random  walk  in  a  plane, 
the  distribution  after  two  steps  shows  singularities,  see  Watson  Sec.  IJ.^S; 
Sec.  15.^6,  eq.  (3);  for  random  walk  on  a  sphere  similar  results  hold,  see 
Roberts  and  Ursell,  1959,  p.  320,  2nd  footnote.) 

For  a  normally  distributed  velocity  field,  a  comparison  of  (2.37)  in 
the  form 


G(k,t)  =  -^   2^  i ,  (2.41) 

(2«)^    o     n:  (n+1): 


with  the  exact  solution  (2.21)  in  the  form 

-    (-l)^(kv^t)2- 

G(k,t)=— ^  2j  TT^ .  (2.i^2) 

(2rt)^   o       2^   nl 


shows  that,  for  short  times^  the  direct  Interaction  approximation  gives  the 
second  moments  of  G(x, t)  correctly.   The  fractional  errors  in  the  fourth, 
sixth,  and  2n   moments  are,  respectively,  1/3,  2/3,  and  [l  -  2  /(n  +  l)  i"]  . 
This  exhibits  in  another  way  the  consequences  of  the  effective  cut-off  in 
"propagation  speed". 


2.5  Solution  for  large  times 

Fot  t  »  t/v  ,  the  diffusing  particles  will  have  suffered  many  displace- 
ments (statistically  almost  unrelated)  from  the  energy  containing  eddies,  and 
we  may  expect  G(x, t)  to  become  close  to  a  Gaussian  distribution.   This  may  be 
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established  directly  (cf.  Taylor^  1921)  by  expressing  the  moments  of  the 
distribution  as  integrals  over  Lagrangian  velocity  correlations.   It  can 
then  be  shown  that,  for  times  large  compared  to  the  Lagrangian  correla- 
tion times  (assuming  these  to  be  finite),  these  moments  approach  Gaussian 
values  corresponding  to  an  eddy  diffusivity 

00 

K..  =     I      df  <u.(f)u.(t  )\^  .  {2M) 


In   a  like  manner,    for  t   »  l/v   ,    the   distribution   satisfying   (2.26)   becomes 
Gaussian  vith  an  eddy  dlffvisivity 


/<.  .   =      /     dt'    /   dx'G(x',t')U.  .(x',t').  (2.i|4) 

ij         J  J      h.       ».  ij  •• 

o 

One  vay  of  proving  this  is  to  extend  to  the  higher  moments  the  arguments 
that  led  to  (2.27).   For  example,  it  is  easy  to  show  from  (2.26)  that 

<^i^jVt>  =  <^i^?  <\^t^  ^  <\\^  <^j^>  ^  <^i^^><^j\> 

-   t  — 

/      dt'(t-t')    /  dx'x.'x'.U,  ,  (x',t' )G(x',t')+5   similar  terms 

J  J  1     J     K-C     ...  ^ 


+  2 


(2.1^5) 
Asstiming  that   expressions   such  as 


i     /     dt'(t-t')    /    dx'x.'x'.U,  .(x',t')G(x',t*) 

tj  J—IJK-O.*,  w 
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converge   as   t  — >  <»,   we   see   that   for  t   »  i/y     the   first  three  terms   on 

2 
the   right-hand  side   of   (2.45)    are   of  order  (iv  t)      while  the   remainder 

are  of  order   t  v  t   and  therefore   are   comparatively  negligible.      Thus  the 
fourth  order  cumulants    can  he  neglected  if  t  »  l/^r   .      Similar  arguments 
hold  in   all  orders   and  show  that  G  approaches  normality.      This   result  may 
also  be   established  by  the   following   alternative  method. 

If  t   »  t/v  ,    the  mean  distance  \/\x     /  the  particles  will  have   tra- 
velled from  their  source    and  will  be   large    compared  to   I   so  that   the 
length  and  time   scales    of  G(x,t)   will  be   large    compared  to   I   and   t/v   , 

M*  O 

respectively.   Under  these  circumstances,  the  only  regions  of  integration 
in  (2.26)  for  which  the  integrand  is  appreciable  are  those  for  which 
G(x-x',t-t')  is  approximately  equal  to  G(x,t).   Thus,  with  K.  .   given  by 
{2.kk) ,   we  have 

2 

SG(x,t)       5  G(x,t) 

_^ =K..-^-^ ,  (2M) 

dt       ij   dx.ox. 
1   J 

since  for  such  large  times  it  is  immaterial  whether  the  upper  limit  of 
integration  in  (2.26)  is  t  or  »,   in  the  isotropic  case  (2.46)  becomes 


SG(x,t) 
~5t 


,2 


kV  G(x,t),  (2.47) 


where  K   is  given  by  (2.51)  with  t  replaced  by  ». 

We  may  regard  (2.46)  as  the  first  term  in  a  series  of  approximations 

based  on  expanding  the  term  G(x-x',t-t')  in  the  integrand  of  (2.26)  in 

a  Taylor  series  about  the  point  x, t,  and  we  may  apply  an  a  posteriori 

«w  ^^^^— ^— — ^— ^— 

check  on  the   reasonableness   of    the   approximation  by  verifying  that,    on 
substituting  in  the  second  term  of  this   series   the  value   of  G(x,t)    derived 
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from  (2.46),  this  second  term  is  small  compared  to  either  side  of  (2.U6). 
We  will  refrain  from  giving  the  analysis,  vhich  is  straightforward.   It 
confirms  that, if  the  relevant  integrals  converge,  the  necessary  conditions 
for  the  validity  of  (2.^+6)  are 

t  »  t/v    ,    X  «  V  t.  (2.48) 

'  o  o 

The   second  condition  arises   from  the   flnlteness   of  the  maximum  propagation 
velocity  which,    as  ve  have   already  seen,    gives   an  artificial   cut-off  at 
the   distance   2v  t.      We  may  expect  that   for  the    actual   case   in  which  this 
cut-off  is  not  present   the   second  of  the   conditions    (2o48)   would  be  unne- 
cessary and  that  the  distribution  of  particles  would  be  Gaussian   at   all 
distances.      In  any  event,   when  the   first   condition  is   satisfied,    it  is 
clear  that   the    fraction  of  particles   affected  by  the   second   condition  is 
negligibly  small. 

It  seems  extremely  likely,    from  the  exact  short-time   result  of 
Sec.    2.2  and  our  present   results   for  large   times,    that   G(x, t)    is  nearly 
Gaussian  for  all  times,    and  that   the  variable   diffusion   coefficients 
defined  in  Sec.    2.5   (cf.    eqs.    2.27,    2.28,    2.31)   will  give  useful  estimates 
of  the  variance   of  the   distribution   at   all  times.      On   assuming  a  Gaussian 
form  for  G(x,t),    we   obtain   from   (2.28)    an  integral  equation   for  K      (t). 
The   approximation  of   (2.35)   by   (2.32)    has  been  proposed  independently  by 
P.G.    Saffman   (unpublished)    in  the  homogeneous    case.      By  assuming   a  Gaussian 
form  for  G(x, t),    and  the  Isotropic  form 
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U.  .(k,t)  = 


k.k. 
1  J 


ij 


exp 


k 
'-  o 


1   P  2^2 


(2.49) 


he  has  obtained  from  the  integral  equation  for  K-  .  the  estimate 

K^0.7v^/k^,   (t-oo),  (2.50) 

vhere  jt   '    K.  is  the  longitudinal  integral  scale  (cf.  Batchelor,  1955^  P'^?)' 

3.    RELATIVE  DIFFUSION 

5-1  Formulation  of  problem;  exact  results  for  short  times 

In  this   section,    we   study  the   correlation  between  the  motion  of  two 
marked  particles  which  are  initially  separated  by  a  distance   small   compaired 
to   I.      The    choice   of  method  is  essentially  that   of  Sec.    2.1   and   again  we 
will  adopt   a  formulation  of  terms   of  Eulerian  moments.     We  introduce  the 
passive   scalar  field  \lf    (x,t)    for  the    first  particle   and  ij/p(y,s)    for  the 
second  particle,    and  we   require  that  both  fields   satisfy   (2.2).      For 
\|f   (x,t),   we   take   as   initial   condition 

1     m. 


^I'/x^t    )   =  5(x-x   ), 

X    iji»       O  «►    ii«0 


(3.1) 


and,    for  ijr   (x,t),   we  take 


tp(y^s    )    =   6(y-y    ) 


(5.2) 
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Let  us  define  a  two  particle  Green's  function  by 

G(x,t;y,s|x  ,t  ;y  ,s  )  =  /i!jx,t)i Jy,s)y  .  (5.5) 

Clearly,  since  the  two  particles  have  identical  properties, 
G(x,t;y,s|x  ,t  :y  ,s  )  =  G(y, s;x,t |y  ,s  ;x  ,t  ) . 

w  „  «.00rt,>.00  a>v>  ~  „UUb/UU 

In  homogeneous   steady  turbulence,    it  depends  only  on  difference  times 
and  difference    coordinates.      We   shall  then  write  it   as 

G(x-x   ,t-t    :y-y  ,s-s    |r   ,T    )   =  G(y-y   ,s-s    :x-x   ,t-t    |-r   ,-t    ) 

«n    wO  O    M  «»o  0    i»0       O  w  «,  O  O    «•   w  o  o      »o  o 

where   r     =  y  -x     and  t      =  s    -t    . 
„o       »o*o  o  00 

It  is   evident   that   if  one   integrates    (3.5)    over  all  x  or  all  y 
the   one-point  Green's   function  of  Sec.    2  is   recovered.      Of  more   interest 
is   the   integral 

^^Z'^'^^^o'lo'K'%^   =/^^^^'*'ri'Hi^o'V2o+5o'%)-  (5.M 

Rfr.t.tlx   ,r   ,t   ,t    )   is  Richardson's   "distance  neighbour  function" 
^w,-'    '    '«o  *o'    o^    o 

(Richardson^    1926).      It   denotes   the  p.d.f.    at  time   t   of  separation  r  for 

a  pair  of  particles  which  at   time   t     were   situated   at   x      and  x     +  r    . 
^  ^  o  wO  <»o       «*o 

For  the  homogeneous   and  stationary  flows  with  which  we  will  be  primarily 

concerned,    R  depends   only  on  r   , t    ,r-r   ,t-t   ,    and  s-s      and  will  be  written 
'  ^  ■^  ,fO     o  -.  »o  o  o 

^^ri:o'*-v^-"oiio''^o)  -  f  ^'^^i-^o'^-^o'^^-^o^^^i-^o^''-%ko''o^- 

(5.5) 
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Let  us  nov  assime  that  the  Reynolds  number  of  the  flov  is  s\ii'ficiently 
high  that  an  inertlal  range  of  vave  numbers,  or  eddy  sizes,  exists.   By 
this  we  mean  that  the  wave  nimibers  which  contain  most  of  the  energy  are  dis- 
tinct from  the  higher  wave  numbers  which  are  responsible  for  most  of  the 
energy  dissipation.   Suppose  now,  that  r  lies  with  this  inertial  range  of 
eddy  sizes.   The  eddies  of  dimension  large  compared  to  r  move  the  two 
marked  particles  together  bodily  without  substantially  altering  the  magni- 
tude or  direction  of  r  .   In  a  frame  of  reference  moving  with  these  large 

^o 

scale  motions,  the  eddies  of  dimension  small  compared  to  r  are  associated 
with  a  small  r.m.s.  velocity  and  have  little  effect  upon  r  .   The  rate  of 
separation  of  the  particles,  in  this  case,  is  dominated  by  eddies  of  dimen- 
sion ~  r  ,  becaiise  such  eddies  make  the  principal  contribution  to  the  rela- 
tive velocity  of  two  points  separated  by  a  distance  r   (cf.  Batchelor, 
1953;>  Ch.  6).   These  eddies  disperse  the  particles  substantially  in  a  time 
of  order 


T(r  )  =  r 
^  o     o 


-  2  n-l/2 

V   -  U. .(r  ,0) 

o       ll^^O 


«lh.  (5.6) 


o 


For  times  short  compared  to  T(r  ),  we  may  apply  argiments  similar  to  those 
of  Sec.  2.   These  show  that  (cf.  eq.  2.9) 


G(x,t;y,s|x^,t^;y^,s^)  =  --^ p(  ^  ,  ^  j  ,       (5.7) 

*"        (t-t^)'^(s-s^)-^   \t-t    s-s   ' 


where  p[u  (x  ,t  ),u  (y  ,s  )lis  the  joint  p.d.f.  for  velocity  u,  at  position 
x^  and  time  t   and  velocity  u^  at  position  y  and  time  s  ,   It  follows 

mO  O  m,  d  ^O  O 
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(cf.    eq.    2.17)    that 


<A  x^  A  y    >   = 


Ax     Ay     G(x,t;y,s|x    ,t    :y   ,s Jixdy, 

1  J  "n  «»  ~0        0   i^O        O        14M    •» 


=   2U.  .(x   ,t   ;y  ,s    )(t-t    )(s-s    ), 


vhere  2^  =  x-x     and  /^  =   j-y   .      Also^    for  these   short  times,    (3'7)    shovs 
that  Richardson's   function  is 


R(r,t,t|x     r     t     t    )    =— i-3    ^(  ^        , 


(3.9) 


where   t/^pvfr  ,t  )"]  is  the  p.d.f.  of  relative  velocity  V  hetveen  the  two 
«M  *.o   o  _j  <** 

points  X  and  x  +  r 

w,0      U.0    U.C 

lows  from  (3.9)  that 


points  X  and  x  +  r  at  time  t  ,  (cf.  Batchelor,  1952a,  Sec.  3)'   It  fol- 

MM  O      •*•  O    u»0  O 


/  A  r.  A  r.  \  =  <V.V.  >(t-t  )^ , 
which,  for  an  isotropic  field  with  r,-axis  along  r  gives 

'  ^  wO* 


(3.10) 


<(Ar^)2>=  <;(Ar2)^>=  2v^^[  l-g(r^,  t^)]  (t-t^)^ 


<(Ar  )2>  =  2v/[l-f(r^,t^)](t-t^)^ 


-\ 


(3»11) 


where  f(r,t)  and  g(r,t)  =  f(r,t)  +  o  ^  Sf(r,t)/Sr  are  respectively  the 
longitudinal  and  transverse  velocity  correlations  at  time  t  for  points 
separated  by  a  distance  r.   Thus,  Richardson's  function  is  initially 
oblate  spheroidal  with  the  line  joining  the  origin  (r  =  O)  to  r  as  axis. 

te-  M>0 
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In  the  particular  case  In  which  P(u^.u_)    is   stationary,    isotropic 
and  Gaussianly  distributed  at  x     and  y  ,    the  double-Fourier  transform 

lanO  uj,0 


p  p      -i(k-x+t-y) 
G(k,t;t,s|r   .T    )   =  -^^       //      e  G(x,t;y,s|r   ,t Jdxdy 


(2n)^ 


(3.12) 


takes   the   form 


G(k,tjt,slr^,T^)    =     ^^^12 


»,       «.      '«.o     o 


(2n)^ 


o   o   o         o   o    o 
v^Tc  t  +v,    4  s  +2n    ,(r  ,t    )k.t.ts 
1  1  ij  «o      o     1   J 


It  follows  that  Richardson's   function 


(3.15) 


ik-r 


R(r,t,t|r  ,0)   =   (2n)^   /   G(-k,t;k,t|r  ,0)6"      dk 


(5.1M 


then  takes   the   form 


R(r,t,t|r  ,0)   = 


«5/2(2,^t)5[l.g(r^)][l.f(r^)jV2 


X     exp 


/      2  2, 

(r^     +  r^    ) 


(2v^t)^    LCl-   g(rj]         [1-   f(r^)]J 


,        (5.15) 


where,    again,    the   r-,-axis  is   along  r   . 


2l|. 


5.2     Intregral  equations    for  G(x,tjy,s|x   ,t   ;y   ,s    ) 
"^ ' ' ' '      '  li^        tfk,  utjO       o  ^*o       o 

When  t-t      and  s-s      are  not   small   compared  to  T(r    ),    it  is   no  longer 

legitimate  to   ignore  the   spatial  and  temporal  variations   of  u.(x,t).      It 

1     m 

is  nevertheless  possible  to  effect  a  partial  summation  of  the  formal  solu- 
tion for  G  containing  terms  of  all  orders  in  the  expansion.   Since  the 
determination  of  G  is  a  problem  which  is  essentially  inhomogeneous  (even 
if  the  velocity  field  is  homogeneous)^  the  Fourier  modes  are  not  weakly 
dependent  (see  Appendix  A),  and  the  methods  expounded  by  Kraichnan  (1959) 
are  not  applicable.   However,  Kraichnan  (1960)  has  recently  generalized 
his  methods  to  inhomogeneous  problems,  and  has  given  an  approximate  equa- 
tion of  motion  for  the  covariance  /ij/(x,  t)\if(x' ,  t '  )  )>  of  a  convected  passive 
scalar  field.   A  straightforward  generalization  to  our  case  of  two  scalar 
fields  yields  the  result 


^  ^^^''^'J'^l^r.'^r.'yr.'^J 


wO   O  ».0   O 


:  /   dt'  /  dx;'U.  .(x,t;x',t' 


SG(x,tlx',t')   aG(x',f;y,s|x  ,t  ;y  ,s  ) 

•►'«••  <*i  wi/i«iOO>tiiOO 


hx.  hx' 

t  1 

o  J 


+/      ds'    /    dy'U.,(x,t;y',s')G(y,s|y',s') 


a^G(x,t;y',s' |x   ,t    ly    ,s    ) 

ilk.  •»»  m^O         O   H>t,0         O 


5x.  Sy'. 
1      J 


(3.16) 


where   G(x,t|x  ,t    )    is   the  one-point  Green's   function  of  Sec.    2.      An  equation 
similar  to    (3.I6)    holds    for  Sg/Ss.     The   equation  for  R(r,t,s|x  ,t  ,t    ,s    ) 

«m^>OmOOO 

that  can  be  obtained  from  (5.I6)  by  integration  may  also,  in  the  case  of  homo- 
geneous velocity  fields,  be  derived  quite  easily  by  a  device  which  reduces 


25. 


the  problem  to  a  fully  homogeneous  one  in  which  the  Fourier  modes  are 
weakly  dependent.   This  matter  is  dealt  with  briefly  in  Appendix  B.   For 
velocity  fields  which  are  homogeneous  and  stationary,  (5>l6)  may  be  writ- 
ten, after  partial  integration,  in  the  form 


SG(x,t;y,s|r_^,T^)     ^2 


at 


°     °  =  V  °s    /  df  /  dx'U.  .(x',t')G(x',f)G(x-x',t-t';y,s|r  ,-r  ) 

ox.  ox.  J      J        »'       XJ    »*  '     ^UM  '4»».        „        '».oo 


^2 
dx. dy. 


ds'  /  dy'U.  .(r  -x+y-y',T  -t+s-s' )G(y' , s ' )G(x,t;y-y' , s-s' | r  ,t  ) 


(5.17) 


The   equation  for  R(r_,t,s|r   ,%    )    to  vMch  ve   referred   above   is 

iMf  aaO        o 


aR(r,t,slr   ,T    )         •   ;s2 
^ »'    '    Vo     o  o 


at 


Sr.dr . 
1      J 


dt'/   dx'U.  .(x',t' )G(x',t' )R(r-x',t-t', sir   ,t    ) 


■-O 


ds'    /    dy'U. ,(r  +r-y',T  +T-s')G(y',s')R(r-y',t,s-s' |r   ,T    ) 


(5.18) 


On  multiplying  this   equation  by  r.r.    and  integrating  by  parts,   we   find 

-^       tJ 


a 


/  R(r,t,s|r   ,T    )r,r.dr  =   2   /  dt'    /  dx'U.  .(x' ,t '  )g(x' ,  t ' ; 

J  lU.  •OOIJ4W  J  JwlJuAU  «uv 


-   2   /    ds'    /    drU. .(r  +r,T  +s ' -t)R(r, t, s ' Ir   ,t    ; 
J  J       to   ij    «o    •     o  ^J  '«o'    o' 


(3.19) 


26. 


By  differentiating   (5-19)   with  respect  to  s^   we  find 


dtSs 


R(r,t,s|r   ,-1    )r.r.dr 


-2   /    drU.  .(r  +r,T  +T)R(r,t,s|r  ,i    ).      (5-20) 

J         m.    IJ    ••O    "<       O  «.  mO       O 


Let  us  take  t  =  0,  s  =0.   Then,  by  (3.20)  and  (2.32),  we  see  that  the 
o       o 

Lagrangian  correlation  for  relative  velocity  V  between  two  particles  which 

■M 

were  initially  separated  by  a  distance  r  is 

UH  O 


<V.(t)V.(s)>       =   2   /    drU.  .(r,s-t) 


G(r,s-t)-R(r-r  ,t,s|r   ,0) 


(3.21) 


This   result   should  be   compared  to  the   exact   result 


<V.(t)V   (s)>       =   /    /    djcdy    <t   (x,t)^   (y,s)u   (x,t)u   (y,s)> 


+  <C'\'r>i:>^>t)ii   {y,s)u  {x,t)u  {y,s)y 


-   <(  >ir   (x,t)\|/   (y,s)u   (x,t)u   (y,s)> 


«*-  ■-    •»» 


<'\|/o(x,t)\|/   (y,s)u  (x,t)u   (y,s)  > 

'^•W  -'-•w  -l-W".  J»«l.  _ 


where   \|r   (x,t)    =  ijf   (x,t|x   ,0)    and  \|f   (y, s)   =   \|/p(y^s|y   ,0)    denote  the  unaveraged 

Green's   functions.      In  stationary  homogeneous   flows,    this   result  reduces   to 

(3-21)    on  the   assimption  that  the   )jf   and  u  fields   are  uncorrelated   (cf.    Sec. 

aw 

2.5). 


27. 


Also,    by  adding  to    (5.19)    to  the   analagous   equation   for  the   deri- 
vative with  respect  to  s,    and  setting  t  =   s,   ve   find 

It  <Vj>  =  /  ^^(;::'^^^i.v°)^iY^  = 


=  ij-   /  dt'    /  dx'U,  .(xSt')G(x',t') 


2  /   df    /   drU.  .(r  +r,t-t') 


R(r,t,t'lr    ,0)+R(r,f,t|r   ,0) 


(3.22) 


For  t  «  T(r   ) ,   R  is  negligible  except  near  r  =  0,    and  the   right-hand  side   of 

O  MM 

(5.22)   becomes 


U.  .(0,0)    -  U.  .(r  ,0) 


t, 


in   agreement  with   (3.10).      For  t   »  l/v   ,    R  is    appreciable   even   at   for  sepa- 
rations  r  of  order   I.      For  such  large   separations,    the    second  term  on  the 
right-hand  side   of   (3-22)    is   quite  negligible,   and  we   find    (cf.    eqs.    2.27, 
^.28) 


<Vj>    =2<x.x.>     =2<yv    >    =kK 


(3.25) 


This  is  consistent  with  the  intuitive  notion  that,  at  such  large  separations 
from  their  source  and  each  other,  the  particles  will  wander  independently  as 


28. 


in  a  true  Browniaxi  motion.  For  Intermediate  ranges  of  t,  it  appears  that 
no  such  definite  statements  can  be  made.  However  a  reasonable  approxima- 
tion appears  to  be  possible  and  this  is  discussed  in  Sec.  3.4. 

5. 5  Solution  of  equations  for  short  times 

For  t  «  T(r  ),   ve  may  assume  that  the  G  factors  in  the  integrands  of 
(3.17)  are  both  negligible  unless  x' ,j' ,x,y   are  small  compared  with  r  .  •  Then^ 

•»  •«<\  Wi      M,  »i«0 

the  U. .   factors   in  the   integrands  may  be   replaced  by  their  values   for  zero 
x',y',x,y^    and  similarly  for  the  time    arguments.      Therefore,    in  isotropic 

*^         tM^         W^     IM 

flows,    we  have 


^G(x,t;y,s|r    ,T    )  .2  ^  ^ 


2  ^ 

^ll(^o^'^o)    ^A^  r    ds'    rdy'G(y',s')G(x,t;y-y',s-s'lr  ,T    ).  (5-2'^) 

±J    mO        O       OX.  Oy  J  J  mu  >^  •«mM<i  hkOO 

•^  o 


Take   a  combined  Laplace   and  Fourier  transform  defined  by 

00  00 

p  p  p        ik-x    p        iJi'y 

G*(k,p;t,q|r^,T^)    =  7-     /   dte'^  /   dse''^^    /    dxe   "  *"  /    dye   *  "'G(x,t;y,s  |r   ,t    ) 

^  o  o 

(3.25) 

Then,   by  (3-2^), 


29. 


pG*(k,p;^q|r  ,T  ) L^G*(t,q)  =  -  ( Sn )  Vv^G*(k,p)G*(k,p;  ^,  q|  r  ,  T  ) 


-  (2n)V  .(r  ,T  )k  t  G*(e,q)G*(k,p;e,qlr  ,T  ), 


(5.26) 


where  G*(k,p)  is  the  Fourier  transform  with  respect  to  x  and  the  Laplace  trans- 
form  with  respect  to  t  of  G(x,t).   By  {2."^) ,    we  have 


(5.27) 


p  +  (2rt)\2v^G*(k,p)  =  ^ 

^    -      (2n)^G*(k,p) 


By   (5.26)    and   (5-27),    it   follows   that 

G*(k,p;t,q|r^,T^)   =  G*(k,p)G*(i,q)/|l  +  ( 2n) V    (r    ,t    )k  t   G*(k,p)G*(^, sj]  . 

(5.28) 

If  we  expand  the   denominator  of  (5.28),   we  obtain 


G*(k,p;i,q|r   ,T    )   =     2-J     ("^^ 


m 


•<M       •>      '«.o     o 


Now,    by   (5-27) 


— im 


(2n)   U.  .(r  .T    )k.t. 
ij'^o'    o'    1    J 


-im+l 


G*(k,p)G*(e,q) 


(5.29) 


G*(k,p)    =   1/  An 


y|^«^[p  +  y 


2  ,,22 

p  +    \/    P     +  '^-k  V 


(5.50) 
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The  inverse  Laplace  transformation  of    [G*(k;p)J  is  therefore   (see,    for 

example^    Watson,    19^^,    Sec.    15-2^    eq.    7) 

(m+1)   J^^^   (2kv^t) 
(2«) 5=^+5   ^^,^f+\   ' 


Thus 


G(k,t;t,s|r   ,T    ) 


WW 

■^     I]     (-lAm+l)  = 


(2n)^ 


k.  t  .U_(r  ,T    )" 
1   n   Ij^    o'    o' 


ktv. 


m 


J     ,(2kv  t)      J        f2tv.s) 
m+l^        1  m+1^        1 


kv  t 


tv  s 


(3.31) 


The  second  moments  of  (3-31)  agree  vith  (3-8);  for  a  normally  distributed 
velocity  distribution,  the  higher  moments  are  given  with  progressively  less 
accuracy  (cf.  Sec.  2.4).   In  Appendix  B,  the  significance  of  (3-31)  in  terms 
of  diagrams  is  discussed.   The  relationship 


f  /,   2^2   2\m-l/2 

^   ,  I4v  t  -  X  ; 
2jt  ml   v^  1       '' 


(lcv^t)^l  ^ 


dx  = 


^V 


,2m+2 


,  if  X  <  2v  t 


,  if  X  >  2v  t 


can  be  obtained  by  induction  from  results  of  Watson  (19^4,  Sec.  13.^2,  eq.  h, 
^  =  1,  using  also  Sec.  13-24,  eq.  l)  or  directly  (Watson,  Sec.  13-l4).   This 
enables  us  to  invert  the  function  G(k,t; t,s I r  ,t  )  in  the  form  (cf.  eq.  2.36) 


51. 


(2n)''(v^t)5(v^s)5     — 


S  (■ 


■1) 


m 


QnH-Dl]' 

[(au)ij' 


Ij'-o'    o'dx  dy. 


G(x,t;y,s|r   ,T   )   =  ( 

«w  IKK)  ••'-'  " 


V 


X 


1  - 


w/t^ 


1  - 


),      2  2 
4v       S 


ni-l/2 


,      if  X  <  2v  t   and  y  <  2v  s. 


m 


if  X  >  2v  t  or  y  >  2v  s.  (3.52) 


As   in  Sec.    2.5,    the  effect   of  the   artificial   cut-off  2v        in  the   p.d.f.    of 
velocity  is    apparent.      It   is    also   clear  that,    if  U      (r   ,x    )    is   small   compared 


to  V  ,  (5-52)  reduces  to  its  first  term 


G(x,t;y,s|r  ,T  )  =  G(x,  t)G(y,  s), 

M    an    mO       O        •»       av 


where  the  right-hand  side  is   given  by  (2.56).      If  the  initial  separation  of 

2 
the  points   is   so   close   in  space   and  time  that  we   can  write  U. .(r   ,t    )    =  v.    6. ., 
■^  ^  ij^»o     o  1      ij 

equations   (5.1^)   and  (5. 51)    show  (cf.   Watson  Sec.    11.^1,    eq.    (12))    that 


R(k,t,sIr^,T    )   =   (2rt)^(-k,t;k,s|r   ,T    ] 

Ml  «»0        O  KM  «»v,         «»0        O 


^         J^[2kv^(s-t)] 


(2«)- 


kv   (s-t) 


or,    inverting. 


32. 


R(r,t,s|r    ,T    )   — >  G(r  +r,s-t);    r     — >  0,    t     — >  0.  (3.33) 


This    shovs.    as  we   expect,    that   in  the   limit   r     — >  0,    t     — >  0,    the   particles 
are  not   separated  by  the   flow  and  simply  move   together   as   one  particle.      (As 
a  further   consistency   checK^    we  note   that   the  Lagrangian   correlation   (3.21) 
is   given   correctly  by   (3-33)    in  this    case.) 

3.^      Solutions   for  large   and  intermediate  times 

For  very  long  times,    t   »  I/y   ,    there   is   a  high  probability  that   the 
particles   have   separated  by  a  distance    comparable  to,    or  greater  than,    i   and 
therefore  wander  independently.      In  fact,    the   solution  of   (3-1?)    is   exactly 
analagous   to  that   of  Sec.    2.5   for  the  one-point  Green's   function,    and  we   find 


G(x,t;y,s|r  ,T    )   — >G(x,t)G(y,s),      (r»t),  (3-3^) 


where   G(x,t)    and  G(y, s)    are   given,    in  homogeneous   flows,    by  the   appropriate 
solutions   of   (2.^4-5).      In  this    case  we   find  R(r,t)    satisfies 

2 
SR(r,t)  S  R(r,t) 

^  —Of/-  ** 

5t  ~  ij        dr.dr .      ' 

1      J 

so  that  R(r,t)    assumes    a  Gaussian  form  corresponding  to   a  diffusivity  twice 
that   characteristic  of  the   one  particle   Green's   function   (cf.    eqs.    2.k'^,    3.23) 

Consider  turbulence   at  high  Reynolds   numbers.      The   existence  of  an 
inertial  range  imples   that  the   spectrim  E(k)    obeys 


55. 


E(k)(ik  =  ^  V   ,    if  k  »  k  , 
^    '  2  o  o' 


E(k)k^dk;  =  |  e/v  ,   if-  k  «  k^, 


where  k  =  (l/l)    and  k  are  vave  numbers  characteristic  of  the  energy  con- 
taining range  and  the  dissipation  range,  respectively,   (e  =  rate  of  dissi- 
pation per  unit  mass;  v  =  kinematical  viscosity.)   It  follows  that 


kE(k)  — >  0,   k  — >oo;   k  E(k)  — >  0,   k  — >  0. 


For  the  later  developments  (cf.  eqs.  5.55^  5' 5^  below),  we  will  require  the 
more  stringent  conditions 


kE(k)  — >  0,   k  — >oo;   k  E(k)  — >  0,   k  — >  0.  (5-55) 


If  the   initial  separation  satisfies 


i   »  r     »  l/k^,  (5.56) 

o  a 


we  expect  that,    in   a  time   large   compared  to  T(r   )   but   small   compared  to 

2  2 

t/v  ,    ^  (r-r   )     y  will  become   large   compared  to  r       but   remain  small   compared 

2 
to   t    .      For  these   "intermediate  times"      (as  we   shall  term  them),    neither  the 

short-time  solution   of  Sec.    5-5  nor  the   long-time  solution  above  Is  valid. 


5^. 


Consider  R(r-r  ,t,t+T|r  ,0)  for  t  >  0.   This  quantity  is  the  p.d.f.  of 

>M  wO         mO 

the  separation  of  two  particles  (released  at  t  =  0  at  a  separation  of  r  ) 

M  O 

one  of  which  has  been  carried  by  the  flow  for  a  time  t,    and  the  other  for  a 
time  t+T.   This  process  may  be  visualised  in  two  stages.   During  the  first, 
of  duration  t,    both  particles  are  carried  by  the  flow.   Their  separation  r 
during  this  time  is  essentially  unaffected  by  the  energy-containing  eddies 
which  give  nearly  equal  displacements  to  both  particles.   It  is  governed 
by  the  motions  (relative  to  the  energy  containing  eddies)  of  dimension  ~  r. 
Now,  in  a  frame  moving  with  the  energy-containing  eddies,  the  r.m.s.  velo- 
city associated  with  these  small  scale  motions  is  very  small  compared  to  v  • 
Thus,  the  mean  square  separation  <^  r  ^  at  the  termination  of  the  first  stage 


2 
is  very  small  compared  to  (v  t)  .   During  the  second  stage,  of  duration  x, 

one  of  the  particles  can  be  considered  as  fixed  in  space  while  the  other  is 

carried  by  the  flow  for  a  further  time  t.   Its  motion  during  this  time  is 

dominated  by  the  energy- containing  motions.   During  the  first  stage  the 

relative  diffusion  is  given  by  R(r,t,t|r  ,0).   Dioring  the  second  stage  the 

••       wO 

further  diffusion  of  the  second  particle  should  be  given  by  G(r,T),  since 
the  energy- containing  motions  are  almost  uncorrelated  with  the  small  scale 
motions.   Thus,  we  expect 

R(r,t,t+T|r  ,0)  =   /  R(r',t,tlr  ,0)G(r-r' ,T)dr' .  (3.37) 

«,       mO  J         m  *iO     «»  >«-     *► 

It  is  clear  that  this  result  is  exact  for  t  =  0  or  t  =  0.   Also,  the  change 
in  the  mean  square  separation  during  the  second  stage  will  be  of  the  order 
of  (v  t)  .   Thus,  when  this  is  large  compared  to  the  value  of  <^ r  ^  at  the 
end  of  the  first  stage,  (3-37)  becomes 


35. 


R(r,t,t+T|r  ,0)  =  G(r,T). 


(5.58) 


Consider  now  (j.lS)  and  its  coianterpart  for  Sr/Ss.   These  shov  that 


aR(r,t,t|r  ,0)     .2 

o 


U  (r',f)-U  (r  +r-rSt') 

Ij   •«  IJ  »0   ^  M« 


X  G(r',f) 


R(r-r',t-t',t|r  ,0)+R(r-r' , t,  t-t ' 1 r  ,0) 


(5.59) 


By  the   argvunents   above,    the   integrand  should  be  dominated  by  contributions 
from  small  t'    (because   of  the  behavior  of  R)    and  from  small  r'    (because  of 
the  behavior  of  G) .      We  vri.ll  therefore  replace  the  R  factors   in  the  inte- 
grand by  R(r^t,t|r   ,0).      It   seems   likely  from  the  preceding     discussion 
that  this   approximation  will  lead  to  results  which  are   at   least  qualita- 
tively correct. 
We  now  find 


SR(r,t,t|r   ,0)  .2 

••  •<  o  o 

5t 


Sr.dr 
1      J    "- 


K      (r,t)R(r,t,t|r    ,0) 

IJ       «  IM.  «0 


(5.^0) 


where 


K.  .(r,t)    =   2    /   dt'    /   dr'    U.  .(r',t')-U.  .(r  +r-r',t') 


G(r',f) 


(5.^1) 


For  short  times  t  «  T(r   ),    (3.^1)   becomes 
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K.  .(r,t)  =  2t 


v^  5.,-  U. .(r  ,0) 
1  xj   ij  »o' 


(5.^2) 


vhlch,  by  (j.^O),  is  in  agreement  with  (J-IO).   For  very  large  times 
t  »  ^/v  ,    there  is  a  high  probability  that  the  separation  of  particles  is 
>  I.      For  these  separations,  the  second  U  factor  in  the  integrand  of  (J-^l) 
is  negligible,  and  the  first  factor  gives 


K.  .(r,t)  =  2  K.  ., 


(5.^5) 


in  agreement  with  the  results  derived  earlier  (cf.  3.25). 

To  calculate  the  form  of  K. .(r,t)  for  intermediate  times,  ve  express 

X  J     m 

(3.ij-l)   in  the   form 


K^.(r,t)   =  2(2it)^   /    dt'    /   dk 


1  -   e 


ik-(r+r   ) 

»         w     "O 


U.  .(k,t')G(k,t' 

X  J      l«  IK 


(3.h; 


Assuming  isotropy,   ve   can  write 


U^.(k,t)=    (5^,  -^ 


E(k,t) 


(3.^5) 


where  E(k,  O)    is   the   energy  spectrum.      It   follows   tha^c   if  we    adopt   spherical 
polar   coordinates    and  write 


R(r-r  ,t,t|r   ,0) 

«,     atlO  mO 


=  R(r,9,0,t), 


(3.^0)   has   the   form 


^R  _  2^    S 
^  "      2  ^ 


"\{t,- 


Sr 
57 


+  \  M(r,t; 


sxn 


a  ,/^ .    ^  sr 

!  sxn  0  ^ 


9  ^  ^' 


1        S^R 
sin  0  d0 


(3.^6) 
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vherej  by  (3.^^)  and  (3.h'^) 


K(r,t)  =  i^(2n)^  /  dt'  /  dk 
o      o 


1   sin  kr   cos  kr 
^    (kr)^    (kr)^_ 


E(k,t' )G(k,f  ), 


i^M) 


M(r,t)  =  2(2jt)^  /   dt'  /  dk 


sin  kr   sin  kr 


cos  kr 


kr 


(kr)^    (kr)' 


E(k,t')G(k,t') 


(3^8) 


For  small  kr^  the  quantities  in  the  square  brackets 


in  the 


integrands   of  (3»^7)   and   {j>.kQ)    are  proportional  to   (kr)    .      Also^    in  the 
inertial  range^   ve  may  take   (cf.    eq.    2.21) 


G(k,t)   =  ?  exp 

(2rt)^ 


1  ,2     2,2 

2  k  v^t 


(3.i^9) 


and   (cf.    Kraichnan,    1959,    Sec.    Q.k) 


E(k,t)    =  E(k)    exp 


1  ,2     2^2 

2  k  v^  t 


(3.50) 


Thus^  since  by  our  initial  supposition  (3'35)j  k  E(k)  — >  0,  k  — >  0,  it  fol- 
lows that  the  integrands  of  (3-^7)  and  {^.k-Q)   do  not  tend  to  infinity  as 
rapidly  as  k   ,  k  — >  0.   Consequently,  the  form  of  E(k)  for  small  k(~'k  ) 
does  not  influence  the  values  of  K(r,t)  and  M(r,t)  appreciably.  Hovever, 
these  quantities  do  depend  implicitly  on  the  energy  containing  range  through 
the  forms  (5.^9)  and  (3-50)  for  G(k,t)  and  E(k,t).  We  will  disciiss  this  in 
more  detail  at  the  conclusion  of  this  section. 
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Three  further  approximations  are  clearly  justified.   First  since 
<^  (r-r  )  /  is  large  compared  to  r   for  the  times  under  consideration, 
r  may  be  neglected  in  (j.^T)  and  (j.^S).   Clearly,  also,  the  particles 
have  in  this  time  lost  all  memory  of  the  initial  orientation  of  r,  so  that 
R(r,0,0,t)  is  spherically  symmetrical  and,  by  (5.46),  satisfies 


5R(r,t)   _  J^    ^ 
~5t  ~      2  ^ 


A(r,t)  55^ 


(5.51) 


p  P 

Lastly,    since  <^r    ^   is   small  compared  to   (v  t)    ,    the  exponential  factors   of 
(5.^9)    and  (5-50)    are   small  at  the  upper  limit  of  integration  over  t  in 
(5.47)    and   (5.ij-8),    and  ve  may  therefore  write 


K(r,t)   =  K(r)    =  k 


dt 


dk 


sin  kr        cos  kr 
— ^  +  ?T 


(kr)- 


(kr)' 


E(k)exp 


1  ,2     2, 2 
V  k  V     t 
3         o 


=    2,/^/ 


dk 


V  k 
o 


sin  kr       cos  kr 
^  +  2 

(kr)^  (kr)t 


E(k). 


(5.52) 


Equation   (3-51)   was  proposed  by  Richardson   (1926)    for  the   Intermediate 
times  discussed  here.      To  investigate  further  the  form  of  the  variable   dif- 
fusion coefficient   (3-52),   ve  will  assume  that  in  the   inertial  range  the 
spectrum  is   a  power  law       of  the  form  (cf-   eq_.    (5-35) 


We  may  now  take  the  wave  number  k     characteristic  of  the  dissipation  range 
to  be 

k,   =  R  ^/^5-n)^ 
do  o^ 


(This  footnote  continued  on  p.  39) 
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E(k)   =  pe""^  v^5-5n  ^-n^      (i  <  n  <  2),  (5.55) 

-where   (3   is   a  dimenslonless    constant   of  the  order  of  unity.      With  this   spec- 
trum,   (5.52)    gives   on  integrating  by  parts  and  use   of  results   of  Watson   (l9^^. 
Sees.    3.h,    15.24), 

^^^)    =  -^r-  ^  ^o  ^      J  n+5/2 

o 


p«v/5       r(l  -  jn)         ^_^       ^  n  n      .        .  .        ,, 

— — —  e  V  r     =Xt  ,    (say).  (5.5^) 

2^+1  n     r(|  +  I  n) 


(Continuation  of  footnote  on  page  58.) 


where 

5  h 

k   =    l/t  =  e/v  -^     ,      R     =  V     /ev  =  V  /k   V. 
o         '  '0^0         o  '  00 

The      intermediate   times"    referred  to   above  may  now  be   defined  more  precisely 

\     /  o  o 

Also,  the  circumstcuices  under  which  (5-58)  is  a  good  approximation  are 

V  -x/l   »  (v  t/l)^^^^-''K 
o       ^  o 

Note  that  (5.55)  may  be  written 
E(k)  =  plv  ^(ki)"". 


ho. 

If  we  suppose  that  at  some  time  t   »  T(r  ) 

R(r,t)  =  ;?(r,t^),  (5.55) 

then  at  subsequent  times  t  (  «  l/v   ),    (3.51)  and  (3-5^)  show  that 


-  i(n+l)  "  •  |(2-n)  -  h{2-nf^^(t-t) 

R(r,t)  =i(2-n)r        ^{^.\)J^Ur  e  '^  l|d^^ 


(3.56) 


where  s  =  (l+n)/(2-n)  and 


o  r   i^5-n)  .  i(2-n) 

jP(bt^)  =  /   r       j^(r,t)jj|r2     )dr  .  (5.57) 


For  t  »  t^ ,    (3.56)   reduces   to 

T,/      ^\                             1  1  -r        /(2-n)    Xt  / -,   ^Qs 

R(r.t)    =  4,^/p  „ r—  .  ^-^  e  /^         ^  ,  (3.58) 

(2-n)^+''/^  ""  r(-l-)  -1- 

^-"        (Xt)   2- 

so  that  „ 

<r2>=    (2-n)V2-n       ^      (Xt)^/^-.  (5.59) 

_ 

The   behavior  of  R(r,t)    and    <^  r    ■)•  given  by   (3.58)    and   (3.59)    is    almost 

certainly  independent   of  the   approximation 

R(r-r',t-t',tlr   ,0)   =  R(r,t,t|r    ,0) 

which  led  from  (3.39)  to  (3.'+0).   In  fact,  it  can  be  shown  that  the  more 

accurate  approximation  method  based  on  (3.37)  and  (3-39)  leads  to  results 

whose  dimensional  forms  are  identical  to  (5.58)  snd   (3.59).   However,  it  is 

worth  confirming  a  posteriori  that  the  approximation  we  have  adopted  is  good. 


(This  footnote  continued  on  p.  ^1.) 
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The  behavior  of  \r  /  as  a  fimction  of  t  is  extraordinarily  sensitive 
to  the  value  of  n  assumed.   Two  cases  are  worthy  of  notice: 
n  =  5/2  (Kraichnan,  1959) 


T^^  ^     1/2   5/2/   1/2      /  2.     2k,     2 
K(r)oC  e  '  r-"    /v     '        ,      <r  >  cC  e  t  /v   , 


(3.60) 


n  =  5/5  (Kolmogrov,  19^^-1) 

K(r)aCe2/5p/5/^^  ^   <  r^  >  ^  A V/  .  (5.61) 

Neither  of  these  agree  with  the  form  proposed  by  Richardson  (1926^  see  also 
Batchelor,  1950): 


K(r)oC  e^/5  ^^/5  ^  /i.2^^  ^  ^^  (5.  2) 


That  this   is   so  is  not   surprising.     Kraichnaji's   direct   interaction   approxi- 
mation does   not  give   an  inertial  range   spectrum  which  agrees  with  that  derived 
from  Kolmogorov's   similarity  arguments.      In  the   same  way,   when  this   approxi- 
mation is   applied  to  turbulent  diffusion,    it  does  not  give   a  diffusion   coef- 

(Continuation  of  footnote  on  p.  ^0.) 

For  this  piirpose,  we  expand  R(r-r' ,t-t' ,t  |r  ,0)  in  a  series  about  r'  =  0, 
t'  =  0;  substitute  in  (3'59);  and,  using  (3.58)7  compare  the  second  term 
of  the  series  with  the  first  (cf.  similar  argument  of  Sec.  2.5).  We  find  that 

this  second  term  is  smaller  than  the  first  by  a  factor  of  the  order  of 

o 
(r/v  t)  so  that,  for  t  «  t/v  (at  which  times  <:^  r  ^  is  small  compared  to 

2 
I   ),\fe   may  expect  (3*58)  "to  be  a  good  approximation  to  R(r,t). 
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flcient  vhich  agrees  with  that  derived  by  these  similarity  argimients, 
even  if  Kolmogorov's  spectrmn  is  assumed.   This  is  because  the  dynamics 
of  diffusion  vlth  a  given  velocity  field  differ  on  the  two  theories.   It 
would"  seem  that  the  behavior  of  \  r  ^  as  a  function  of  t  might  provide  a 
sensitive  test  by  which  to  confront  with  experiment  different  assimiptlons 
about  the  structure  of  the  inertlal  range. 

Kralchnan  (1959)  has  made  a  detailed  comparison  between  the  direct- 
interaction  approximation  and  the  Kolmogorov  theory.   He  has  traced  the 
difference  in  the  Inertial  range  spectra  to  the  difference  in  the  role 
played  by  the  energy- containing  eddies  in  the  two  cases.   In  Kolmogorov's 
theory  these  eddies  merely  convect  the  small  scale  motions  without  influ- 
encing their  dynamics^  whereas  in  the  direct-interaction  approximation 
this  is  not  so.   In  the  same  way,  on  arguments  of  the  Kolmogorov  type, 
the  relative  diffusion  of  particles  should  be  independent  of  the  energy- 
containing  eddies.   However,  the  application  of  the  direct-Interaction 
approximation  has  led  to  results  which  depend  on  the  energy- containing 
motions.   If  we  wish  to  modify  our  formalism  in  such  a  way  that  our  results 
depend  only  on  the  small  scale  motions,  we  would  transform  to  a  frame  of 
reference  moving  with  the  energy- containing  eddies.   We  would  expect  that 
results  of  the  form  (J.-^O),  (j-^l)  would  be  qualitatively  correct,  pro- 
vided that  G(k, t)  now  described  one-point  diffusion  relative  to  a  source 

moving  with  the  energy- containing  eddies.   On  similarity  argijments  of  the 

^  2/3  1/3 

Kolmogorov  type,  G(k,t)  would  then  depend  only  upon  k    e    t.   Similarly 

E(k,t)  would  describe  the  structure  of  the  small-scale  eddies  in  a  frame 

of  reference  moving  with  the  large  scale  motions,  and  would  take  the  form 


i^5. 


E(k,t)    =  E(k)f(k^/V/^t). 


1/2 
The   overall  effect  of  these  modifications  vould  be   that    [kE(k)J  would 

appear  in  place   of  v     in  the  expression   (5-52)    for  K(r),     Kraichnan   (1959^ 

Sec.    9.1)   has   shown  that  this   substitution  resolves   the   conflict  between 

the  Kolmogorov  theory  of  turbulence   and  the  theory  based  on  the  direct-inter- 

1/2 
action  approximation.      The  quantity    [kE(k)3  may  be   considered  as  the   r.m.s. 

velocity  associated  with  the  motions   of  wave  numbers  k  as  they  are   convected 

by  the  large-scale  motions.      Substitution  of  this   quantity  for  v^  in   (5«52) 

gives 

|(n-l)       i(5-5n)     |(n+l) 
K(r)=<:e^  v/  r^  ,  (3.63) 

which;,  taking  n  =  5/5^  leads  to  (5.62).  With  these  changes  in  the  inter- 
pretation of  G(k,t)  and  E(k,t)^  the  sensitivity  of  our  results  to  the  form 
of  the  inertial  range  spectriim  remains. 
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APPENDIX  A: 

Derivation  of  equation  (2.26) 

When  the  velocity  field  is  spatially  homogeneous,  the  probleTn  of 
diffusion  fi'ora  a  point  source,  although  apparently  possessing  only  radial 
symmetry  even  in  the  isotropic  case,  can  always  be  rephrased  as  a  homo- 
geneous problem.   For,  since  the  equation  (2.2)  is  linear,  the  response 
of  the  system  to  an  initial  disturbance 

^Kx,t^^)  =  e^^'-   ,  (A.l) 

is 

<t(x,t)>  =Jdx^G(x,t|x^,t^)e^^-5o   ^  (A. 2-) 

and.  since  for  a  homogeneous  velocity  field,  G(x,t|x  ,t  )  depends  on  x  and 


x   in  the  combination  x  -  x  only,  equation  (A. 2)  -can  be  rewritten 
WQ  ~    -^o 

<^(x,t)>  =  (2«)5  G(k,t|t^)e^5-25   ,  (A. 3) 


where 


^  r  -ik.(x-x  ) 

'G(k.t|t  )  =  — ~      /  G(x-x  ,t|t  )e  ""  "  "^   d(x-x  )         (aA) 
-^         o    ^2«)5  J    -  ~o   I  o 


Lon 


is'the  Fourier  transform  of  the  Green's  function  G(x,t|x  ,t  ).   Equatic 
(a. 3)  proves  the  average  response  matrix  of  the  Foiu-ier  modes  is  diagonal 
(when  the  velocity  field  is  homogeneous)  and  that  {2n)      G(k,t|t  )  is  the 
average  response  funct.Lon  for  mode  k. 

Having  established  this  correspondence,  we  will  now  derive  an  approx- 
imate equation  for  the  response  fiuiction  by  a  method  parallel  to  that 


^5. 

employed  by  Kraichnan  (1959)  ^ov   the  velocity  field  response  function. 
For  simplicity^  we  will  suppose  henceforth  that  the  velocity  field  is  also 
statistically  stationary.   The  modifications  necessary  if  the  field  is  not 
statioiiary  are  easily  included. 

It  is  convenient  to  Introduce  the  artifice  of  cyclic  boundary  con- 
ditions over  a  large  cube  of  side  L  in  order  to  expand  i|f(x,t)  and  j^(x,t) 
in  Fourier  S'Jims  rather  than  Fourier  integrals: 


^(x,t)  =  Zj     t(k.-t)e^-'-  ,  [j(k)  =  r(-i5T]  ,  (A. 5) 

u(x,t)   =  X!    ^(k.t)5^-i^-^   ,   fuCk)   =  u*(-k)]      ,  (A. 6) 

(cf.    K-^    eq.    2.1)  0      Equation  (2.2)   may  be  written 

M^.tl   ^    .        V^    ll(p,t)t(q,t)    =    0   ,  (A. 7) 

where  p  =  k-q.  The  response  function  g(k,t)  for  mode  k  is  the  solution 
—   -»-■_>  —  " 

of  (a. 7)  under  the  initial  conditions 


l}(k,t)  =  g(k,t)  =1 


t(q.t)  =  ^.  (q,t)=  0  ,  q  7^  k 


.  t  =  0  ,  (A. 8) 


This  paper  will  be  designated  by  'K.'  in  these  Appendices.   Some  differ- 
ences in  notation  should  ha  noted:   In  K.,  g(k,t)  refers  to  the  velocity 
field  (impulse)  response,  function  and  not  to  the  response  function  for 
(a. 9)  below.  Also,  in  K.,  g  refers  to  an  averaged  response  while,  in  this 
appendix,  it  does  not;  the  average  being  denoted  by  <^g'^.  Further,  the 
notational  distinction  "between  a  quantity  and  its  Fourier  transform  is  dif- 
ferent from  tiiat  adopted  in  this  paper. 


1^6, 

(cf.K.  Sec. 2.1).  By  the  equation  of  motion  tj^(q,t)  and  the  direct  inter- 
action approximation,  we  find  (cf .  K.  eq.  2.2*1)  that  for  solution  (A. 8) 

t 

■^Jq.t)  =  -iq.  f    u  f-p,f)g(k,f)g(q,t-t')dt'   ,  (A. 9) 

t,  -  f-  J         ^      ^  ~      - 

o 

Tnus  by  (2.50)  we  nave 

t 
^%^  =  -  X!  k.k.  r  u.(p,t)u.(-p,f)g(q,t-t')g(k,f)df  ,  (A. 10) 


and,  on  averaging,  using  the  principle  of  weak  statistical  dependence 
(cf .  K.  Sec.  2.2  and  eq.  (2.25)), we  find 

t 

|^.<g(k.t)>=  -i;  k.k  /<u.(-p,t)u.(-p,t.)> 
^        o 

<g(q,t-f)>$(k,t'))dt'  .    (A. 11) 

Now  let  us  take  the  limit.   Make  the  transition  L  — >  go  .   Let 


=  Limit   /  LVyu.(k,t)u.(-k.t')V 
L— >  GO  [2nP   ^  ^    J     -y        / 


U..(k,t-t')  =  Limit   fj^y^C^,iKt)u.{-k^t'jy,  (A.  12) 


so  that 


U.  .(x-x',t-t')  =  /u.(x,t)u.(x',t')\  =  Tu.  .(k,t-t')e^-"^^'-'Vk  ,   (A. 13) 


and  let 

5 

that  equations  {A.h,    5,  8)  are  consistent  (cf.  eqs.  5-2,  3-3) ■     Then 


r(k,t)  =  Limit    (^y<£ih^)y  ^^-^^^ 

L  — ^>  CO  \   / 


so 

t 


5G^^=  -(2«)5  k  k.  /  df  /  dqU..(p,f)G(q,t')%t-f)   .       (A. 15) 


^7. 

This  result  can  be  returned  to  physical  space  by  vriting  it  as 

t 
M|^^^  -k.k.  r  df/dp  /dq  /L-  U..(p,f)G(q,f  )G(k,t-f)e^-  •^£"'3."-^  ; 

(A. 16) 

i.e.  (of.  A.-^,  13) 

t 

^It^  =  -k.k.  r   dt'  Tdx'  U.  .(x',f)G(xSt')G(k,t-tOe"^--^'  •    (A.l?) 
o 

Using  (A.k)   again,  this  can  be  rewritten 

2      ^ 
-'^-'-^  =  ^-  /  ^'^'   /^'  U  (x',f)G(x',f)G(x-x',t-tO  .      (A.18) 

o 

We  now  give  an  alternative  derivation  of  (A. 15)  based  on  su-iming 
selected  terms  from  all  orders  of  the  foruial  expansion.   For  simplicity 
we  will  treat  the  one  dimensional  case  only.   Then  (2.2)  is  replaced  by 

-%^  =  -  ^  [u(x,t)t(x,t)]   .  (A.19) 

On  making  expansions  analagous  to  (A. 5,  6),  we  obtain  from  (A.19) 

^|^=  -ikXl  u(p,t)t(q,t)   ,  (A. 20) 

q 

where  p=k-q.   (We  will  again  suppose  the  mean  flow  a(o,t)  is  zero.)  The 
formal  solution  for  the  response  function  for  mode  k(cf .  A. 8)  is 


^8. 


g(k,t)  =  1  -  y]  k(k-l)[l*u(0*u(m)]. 
t+m=0 


+  ±y]       k(k-t)(k-l-m)[l*u(l)*u(m)*u(n)]^  ^^_2i) 


+       N""    k(k-^)(k-^-m)(k-^-m-xl)  [l^(.u(e)*u(m)*ii(n)-x-u(r)3^ 
t+m+n+r=0 


where  the  *   denotes  the  process  of  convolution:  - 

t 
[f*g]^  =f    f(t')g(t-f)dt'   . 
o 

To  determine  the  response  'if   (q^t)  induced  in  mole  q  by  the  initial  exci- 

k 

tation   in  mode  k,    we   can  write    (A.20)    in  the   form 


\{q,t)    +   IqJ    df    u(-p,f)i(k,t')    =    -i<iyj    /    u(l,t')^^(q-t,f)cif, 
o  t^-p     o 

(A. 22) 

and  iterate  starting  with  a  zero  "  approximation  derived  by  setting  ^jr  =  0 

on  the  rlght^  followed  by  a  first  approximation  obtained  by  substituting  the 

th  "^ 

zero   approximation  for  t,  into  the  right-hand  side  of  (A. 22),  and  so  on. 

k 

In  this  way,  we  find 

u(p,t)tj^(q,t)  =  -iqu(p,t)(l*{u(-p)g{k)}]. 

-qu(p,t)  ^   (q-0[l*u(^)*{u(-p-^")g(k)}]^        ^^23) 


't 


+iqu(p,t)  y~'     (q-t)(q-t-m)[l*u(e)*u(m)*  {u(-p-t-m)g(k)}]. 
t,xn 


^9. 


where  the  suras  exclude  indices  for  which  I,    t  +  m,    t  +  ra  +  n,  are 

equal  to  -p. 

So  far,  the  analysis  is  exact.   Now  we  make  a  selection  of  terms. 

Equation  (A. 7)  shows  that  to  each  interaction  which,  through  the  agency 

of  'u(p),  induces  a  rate  of  change  in  t(k)  proportional  to  i|A(k-p),  there  is 

a  coii.jugate  interaction  which,  through  the  agency  of  u(-p),  induces  a  rate 

of  change  in  iif(k-p)  which  is  proportional  to  »|^(k) .   There  are,  of  course, 

less  direct  couplings  between  i|^(k)  and  ■lir(k-p)  which  involve  more  than  one 

intermediate  u  mode,  but,  in  the  formal  solution  (A. 21)  for  g(k,t)  and  the 

exTDcessioii  (A. 25)  for  u(p,t)  f,  (q.t),  we  will  only  retain  a  term  if  all 

it 

the  interactions  of  which  it  is  composed  occur  in  conjugate  pairs.  (We 
note  that  this  automatically  excludes  all  terms  with  an  odd  n^omber  of  u 
factors. ) 

We  can  also  explain  our  selection  r-ule  by  a  diagramatic  i-epresentation 
of  "che  terms  composing  the  series  (A. 21).  A  typical  term  containing 
2n  u-f actors  will  be  said  to  be  of  order  n.   The  wave  numbers  of  the  \|f(k-p) 
modes  involved  in  the  interaction  with  ^{'k.)   form  the  segments  of  the  case 
line  of  the  diagram  (see  figure  l)  and  the  ingoing  arrow  on  the  left 
marked  k  represents  the  mode  t(k)  itself.   The  next  segment,  marked  (k-i), 
represents  the  \|f(k--t)  mode  with  which  it  interacts,  the  n{l)   mode  concerned 
in  the  interaction  being  represented  by  the  loop  marked  I   leaving  the 
base  line.  The  following  Intersections  of  loop  and  line  have  a  similar 
meaning.  The  selection  rule  ensures  that  to  every  vertex  A  involving 
i|f(s),  u(l),  'ijf(s-i)  there  is  a  conjugate  vertex  a'  involving  ^{s-l),   u(-t), 
i|f(s)  .   Thus  the  sum  of  all  the  u-factors  involved  at  vertices  between  A 
and  A  is  zero;  i.e.  there  is  an  even  number  of  such  vertices  and  these 
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m 


.C^■^■ 


k  k-l 


k  k-m  k 
(a) 


k  k-l 


k-l-m  k-m  k 
(c) 


Fig.  1:   Diagrams  for  n  =  2  representing  three  non- 
zero terms  of  [l*u(  e)*u(m)*ii(n)*u(r)]   of 
equation  (A. 21).  (a)  t+m=n+r=o, 
(b)  ^  +  r  =  m  +  n  =  0,  (c)  i  +  n  =  m  +  r  =  0. 
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are  conjugate  in  pairs.   Thus  the  linss  joining  conjugate  vertices  do 
not  cut  each  other. 

Two  consequences  of  the  selection  rule  should  be  noted.   First,  the 
only  terms  of  the  formal  expansion  which  are  Included  are  those  which 
involve  pairs  u(^),  u(-t).  These  terms  are  the  same  in  normal  and  non-nor- 
mal distributions  having  the  same  covariance.   Secondly,  not  all  terms 
involving  pairs  u(t),  u(-l)  are  included;  for  they  will  not  be  conjugate 
vertices  if  the  loop  joining  them  cuts  the  loop  joining  another  pair. 
For  example  in  figure  1  all  the  Interacoions  are  shown  for  n=2  and,  by 
the  selection  rule,  we  retain  only  l(a)  and  l(b). 

Consider  now  a  general  term  in  the  series  (A. 23). 

iq(q-^)(q--^-m) u(p,t)  [l^u(t)*u(m)*u(n)x- *{u(-p-e-m-n- )g(k)}J  ,  . 

The  initial  u(p,t)  terra  must  find  a  mate  u(-p)  within  the  convolution  pro- 
duct.  Suppose  it  is  not  u(-p-i-m ).   It  cannot  be  any  of  the  even 

placed  terms  u(m) ,  for,  if  it  were,  there  vould  "be  an  odd  num- 
ber of  vertices  on  the  diagram  between  ■u(p)  and  u(-p)  and  these  could  not 
all  find  mates  without  cutting  the  loop  joining  u(p)  and  u(-p).   The  mate 

for  u(p)  must  therefore  be  one  of  the  odd  placed  terms  u(t),  u(n), 

There  is  then  an  even  number  of  vertices  on  the  diagram  between  u(p)  and 
u(-p)  and  these  vertices  can  and,  by  the  selection  rule,  must  be  joined 
without  cutting  the  loop  joining  u(p)  and  u(-p);  i.e.  the  su>n  of  the  wave 

numbers  I   +m  +  for  vertices  lying  between  u(p)  and  u(-p)  must 

be  zero;  in  other  words,  the  sum  of  these  Intermediate  wave  numbers  together 
with  the  wave  number  -p  for  the  mate  of  p  itself  must  equal  -p.   But  these 
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terms  are  pirecisely  the  ones  excluded  in  the  s'jjiniiatlons  of  (A. 26).   Hence 
we  have  reached  a  contradiction.   Thus,  the  mate  for  u(p)  must  be 

u(-p-t-m.  .  . )  and  the  sum  ^+m+ of  all  Intennediate  wave-numbers 

must  be  zero.   Now  compare  the  general  term  of  the  series  (A. 25)  for 
u(p,t)  -ijf  (q_,t)  with  the  corresponding  term  [as  derived  from  (A. 21)]]  for 
the  series  for 

-iq  u(p,t)  Q  u(-p)g(k)}  *g(q)]^ 

It  will  be  readily  seen  from  the  above  argument  that,  under  the  selection 
rule,  the  two  series  are  identical;  i.e. 

t 
u(p,t)tj^(q,t)  =  -iq  r  u(p,t)u(-p,t')g(k,t')g(q,t-t')dt'. 

o 

Hence,  by  (A.20),  and  the  principle  of  weak  statistical  dependence  (K.  Sec.  2.2), 

t 
a<gCp}>^_^-^  l/  'it'<u(p,t)u(-p,f)><g(k,t')Xg(q.t-f)>,   (A. 21^) 
q     o 

as  before.   (Tnis  proof  generalizes  easily  to  5-dimensioas . ) 

As  further  confirmation,  we  can  show,  by  the  direct  counting  of  dia- 
grams, that  for  large  k(»  i        )  the  selection  criterion  gives 

<;g(k,t)>  =.  J^(2kv^t)/kv^t   ,  (A. 25) 

which  (cf.  eq.  2.57)  is  also  the  solution  to  (A. 24)  under  these  circum- 
stances.  For  large  k,  all  the  terms  k--E-,  'k-l-m,    etc.  of  (A. 21)  are  equal 
to  k  with  negligible  error.   Also,  any  convolution  term  of  order  n  is  asymp- 
totically equal  to 
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(v^t)2"  /  (2n): 


(since,  as  we  can  shov  a  posteriori,  <!^g(k_,t)^  is  negligibly  small  unless 
t  «  t/v  )  .   Hence 


o 

CO 


2n 


<g(k,t)>     =2     (-1)"  N^>Y^t)^7(2n):  (A. 26) 


where  N     is  the  mamber  of  non-intersecting  diagrams   of  order  n. 
n  CD 

To  determine  N  ,   we   introduce  the   idea  of  a  block.     A  diagram  will 
be   said  to  be   composed  of   s  blocks   if  there  exist   s  particular  conjugate 
pairs   A,  A      ,   A  A.,      A^A        with    riie  following  property:    that  both 

-L   J-       £_  1_  s   s 

I  I 

vertices  of  every  other  conjugate  pair  B  B  lie  between  A  A   for  some 
choice  of  r.   In  other  words,  there  are  no  loops  in  the  diagram  which 

r 

(toj^XDloglcally)  lie  above  any  of  the  pairs  A  A    .   For  examples,  see  fig- 
ure 2,  in  which  the  s  loops  composing  the  outline  o?  the  blocks  are  identi- 
fied by  heavy  lines.   Suppose  the  number  of  diagrams  ol"  order  n  composed  of 
s  blocks  is  N    (l  <  s  <  n)  .  We  proceed  by  inductioa.   We  consider  an  s 

Xl^  s 

block  diagram  of  order  n-1,  and  add  an  additional  loop  whose  left-hand 
vertex  L  lies  to  the  left  of  the  entire  s  block  diagram.   There  are  (s+l) 

possible  positions  for  the  right-hand  vertex  R:  P,,P,-, P  -, 

(see  figure  5)>  and,  if  we  choose  P_,  the  diagram  of  oMer  n  thereby  gene- 
rated  has  r  blocks.   Remembering  that  L  is  always  to  the  left  of  the  entire 
diagram  of  order  (n-l),  it  is  not  difficult  to  see  that  every  diagram  of 
order  n  can  be  generated  without  duplication  by  applying  this  process  to 
every  (n-l)  order  diagi'am.   Since  an  (n-l)  order  s  blook  diagram  produces 
one  diagram  each  of  p  blocks  (p=l,  2 ,  s+l),  we  see  that 


5^. 


.z:x 


zx 


n=3,   s=2 


/^^r\r^ 


fr^^r^ 


r\r\ 


n=5,  s=  I 


n  =  7. 

Fig.  2:   Illustrations  of  block  structure  of  diagrams. 
Outlines  of  "blocks  are  identified  by  heavy  lines. 


s  =  3 


(^TA.tfTA     iflTh  ffTTh^        (fTTh  (f77\. 


L  P, 


s+i 


^s-i  Pf+l 


Pr. 


Fig.  5:   The  generation  of  a  diagram  of  order  n  from  a 
diagram  of  order  n  -  1.   Shaded  areas  denote  blocks. 
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n-l 

N    =  V^  N  ,  (A.?7) 

n.s    ^   n-l,r  , 

r=s-l 


with  the  convention  N  ^  =  0  .   It  follows  that 

n,0 


-^r^^H   ,  (A.28) 


n,r    nl(n-r; 


and,  since  every  n  order  diagram  gives  rise  to  one,  and  on.ly  one,  n+1 
order,  1  block  diagram,  we  nave 


n    n+1 , 


Therefore,  by  (A. 26)  , 


>2n 


00        (kv  t,)  -    J  (2kv  t) 
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APPENDIX  B: 

Derivation  of  equation  (3.I8) 

As  in  appendix  A,  we  will  suppose  that  the  velocity  field  is  homo- 
geneous and  stationary^  and  consider  homogeneous  forms  of  t-,  and  ^      .      We 
again  adopt  cyclic  boundary  conditions  and  expansions  of  the  fonn  (A.5j  6); 

ik.(x-x  ) 
t^(x,t)  =2  "^iCk^t)  e °    ,  (B.l) 

k 
k 


u(x,t)   =  y^  u(k,t)  e^--~        ,  (B.3) 

k 


(cf.  K.  eq.  2.1  and  footnote  to  Appendix  A).   Then  equation  (2,2)  for  the 
\|f  field  may  be  written 


S^n(k,t)  -ip.X 

^      +  ik.  y^  u.(p,t)  e  --'-*°  t.(q,t)  =  0   , 


ix  ~o 
at 

<1 


(B.M 


and^  for  the  f  field. 


-— ~-+  ik.  2  ^(P^^)  e  --  ^2 


ik,  V  u  rp,s)  e  ~~°  to(q.s)  =  0  ,  (B.5) 


where  p  =  k-q.   By  (B.^)  we  have  (cf.  K.  eqs .  5-^^  5.5) 


ip.x. 


St 


<Vl(-k,t)  t2(k,s)>=  ik.  ^  ^u.(-p,t)y^(-q,t)  t2(k,s)>e  "^  ""^.  (B.6) 


57. 


Nov,  in  the  limit  L  — ^>  oo ,  the  direct  interaction  approximation 
(K.  Sec.  2.k)   gives 

(B.7) 
where  (cf.  eq.  A. 9  ani  K.  Sec,  j) 


r  ^  -ip.x 

T   ,  (-q,t)  =  /   iq  u.(p,f)t,(-k,f)g(q,t-f)e  "^ -°  dt '   ,   (B.8) 

t 
o 


tp   (k,s)  =  -/   ik  u  (p,s')i(r  (q,s')g(k,s-s')e  -"""ds'   .     (B.9) 
2,0^-'  J  JJ~     ^~      -^ 

•^        s 
o 


Thus,  applying  the  '-/eak  statistical  dependence  principle  (K.  Sec.  2.2), 
we  find  that,  in  the  limit  L  - ->  oo,  (cf .  eqs .  A.  12,  13,  1^4- ) 

^-  R(k,t,s  Ir  ,T  )  = 
<jt  ^^'    '    u o'  o 

t 
-(2n)5  k.k.y^  dfy^dq  U..(p,f)G(q,t')R(k,t-t's|r^,Tj  (^^^^^ 

o 


+(2::)'^  k.k.  /  ds'  /  dq  U.  .  (p,  t  -t+s-s  '  )e  ~""^  a(k.s  '  )R(q.t,  s-s  '  Ir  ,  t  ) 

1  J  J        J         "      2-^    -V   O  ~J'        w  '—"O   o 


Returning  to  x-space  by  the  arguments  of  (A. 15)  to  (A.18)  this  becomes 
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^-  R(r,t,s Ir  ,  T  )  = 


/• 


,2 
J 


t 


^^—  1  /  dt'  /  dr'  U.  .(rSf  G(r',t'  R(r-rSt-t' ,s  r  ,t  )      (B.ll) 
1  ,1   "-^^     ^ 


-        ds'  /  dr'  U.  .(r  +r-r',T  -t+s-s')G(r',s')R(r-rSt,s-s' Ir  ,T  )  I 
I  ,  /   — '   ij  "o  --  ~  '  o       '    ^^  '      '    ^^  ^  '    '  1^0  o  r 

o 

J 

which  is  equivalent  to  (3.l8). 

An  alternative  derivation  of  (B.ll)  can  be  given  which  is  based  on 
Slimming  selected  terms  from  all  orders  of  xhe  formal  expansion  of  R.   We 
will  not  discuss  this  method  in  detail,  but  we  will  describe  how  it  differs 
from  the  similar  method  of  appendix  A,  and  for  this  purpose  we  will  once 
more  discuss  the  one  dimensional  case.   The  formal  solution  of  (B.i^■)  under 
the  initial  condition  \|/  (x^t)  =  5(x-x  )  is  given  by 


L5r_j_(-k,t)  =  1+i  V  k[l*u(t)]^e    °  -  V  k(k+^)[l*u(^)^^u(m)]^e        ° 

I  t,m 

(B.12) 

-i(t+m+n)x 

-i  "V^  k(k+t)(k+t+m)  [l*u(e)*u(m)*u(n)]^e  °+ 

^,m,n 


5  "^  /    \ 

Multiply  this  series  term  by  term  by  the  similar  expansion  for  L  A;   {\,\,) 

and  average,  retaining  only  conjugate  interactions .   Now,  however,  the 

mate  of  'u(p)  in  one  series  may  occur  in  a  term  of  the  other  series,  and 

it  is  precisely  such  interactions  which,  on  the  direct  interaction  approxi- 
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mation,  give  rise  to  the  correlation  between  ^^    and  to-  ^^   ^^  Appendix  A, 
a  diagramatic  representation  may  be  made  of  the  formal  expansions.   Nov^ 
however,  there  are  two  base  lines,  one  for  ^     and  one  for  ^|p  .   The  u  modes 
are  represented  by  loops  Joining  conjugate  vertices  either  on  the  same  or 
on  different  base  lines,  and,  as  in  Appendix  A,  these  loops  cannot  cross. 
A  typical  term  of  the  series  for  \|r  (-k)t2(^)   is  shown  diagramatically 
in  figure  k:      it  is  obtained  from  the  product  of  a  second  order  term  In 
\|f-|(-k;  and  a  fourth  order  term  In  tp(kj  . 

On  the  basis  of  this  selection  rule  of  Appendix  A,  we  will  now 
recover  the  solution  of  Sec.  3-3  for  large  k.  Under  these  circumstances, 
all  the  terms  k  +  I,    k  +   I   +  m,    etc.  of  (B.12)  are  equal  to  k  with  negli- 
gible error.  Also,  a  pair  of  conjugate  vertices  will,  asymptotically, 
give  a  contribution  of 


2^     u(-t,o)u(t,o)  =  v^   , 


-^^  -y^ 


if  both  vertices  refer  to  i)r  or  ^^^   ,  and  give  a  contribution  of 


1]  u(-t,t^)e''"°u(l,s^)e  ''''''  =<u(x^,t^)u(y^,s^)>=U^^(r^,x^), 
t 


^^ 


if  one  vertex  refers  to  t-  and  the  other  to  to  •   Henct 


<;^t^(^k-,t)t2(k,s))>  = 


yp  .  (--) 


m+n 


/   I   (t+2m)l(t+2n): 
n,m,  I 


11^  o,  o 


2 

V, 

1 


^    ,(kv  t)^+2m(kv  s)'-^2n^   (B.lj) 
m;n;t   1         1 


6o. 


2      k      k+l  k+I+n         k+l+n-p  k+l+n-p-q   k+l+n-p   k+l+n 

Fig.  h-:      A  typical  term  in  the  formal  expansion  of  <^\|;  to/* 


Q.+I    Qt      Qm^     ^ 


Fig.  5i   A  general  term  <("  Iji'  ij/^  /  showing  location  of  termini 
of  links.   Blocks  are  shaded. 
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whsre  N    .  is  the  number  of  non-intersecting  diagrams  with  I  +  2ni  points 
m;n;  t 

on  the  t-|  -has-*  line,  t  +  2n   points  on  the  tp  base  line  and  I   limes  between 
them. 

To  determine  N    ,,  we  first  evaluate  the  number  N(iTi,r;  n,  s;  l) 
m^n^  l 

^^ 
of  diagrams  in  which  the  t  +  2ra  points  lying  or^   the  t  base  line  are 

composed  of  r  blocks  and  I   lin}:s,  and  the  t  +  2n  points  on  the  t^  base 

line  are  composed  of  s  blocks  and  I   links.   For  example,  figure  h   shows 

one  diagram  counted  in  N(2,1;  h,    1;2).  There  are  (r+l)  possible  loca- 

tTons  P  ,  P P  ,  for  the  termini  of  the  I   lial<.s  on  the  -ijf,  base 

1'  2.'  r+l  1 

line  and  (s  +  l)  possible  locations  Q^,  Q^,  Q„_^-^^  foT  the  termini 

on  the  "ijAp  base  line  (see  figure  5),  and  any  number  of  termini  can  be 
placed  at  any  of  these  locations.   There  are  therefore  (r  +  t)  .'  /  rl  l\ 
ways  of  placing  the  I   termini  on  the  i|f,  base  line  and  (s  +  l)\    /   si  l\ 
ways  of  placing  the  other  I   termini  on  the  to  base  line.  There  is  only 
one  way  of  joining  the  I   termini  on  each  base  line  with  non- intersecting 
links.   Hence,  in  all,  there  are 

r  ft :     s  :  I ! 

ways  of  linking  the  base  lines.   As  we  have  seen  (cf.  A. 23),  there  are 

r(_an-_r2ll'    s(2n-s-l)  ! 
m:(m-r):'    n:(n-s): 


ways  of  forming  the  blocks  on  base  lines.   Hence. 
N(m,r;  n,s;  I)      =  M(m,r;  l)   M(n,s;  l) , 


where 
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../     ,\    (r+l):   rygrn-r-l).' 
^  -  '  '    r 1 1 1    mi  (m-r j  : 


(B.lM 


It  follovs  that 


N 


m;n;  I 


m 


11 

s=l 


Mi'n^s^  t) 


(B.15) 


Now^  let 


m,  t 


m-l 


^.'m:  V^  M(m,r;l)  =  \  ^ 


_  ^- n  (mft-r):   (m+r -1)  ! 


rl      (m-r-1 


r=l 


From  the  calculations  of  Appendix  A  (cf .  A. 28;,  A. 29),  we  know  that 


(B.16) 


S   ,  =  (2m-l);  /  ml  ,  S   ^=  (2m):  /  (m+l)  !  , 
m, -1   ^     '      '  •"   m,  0 


(B.17) 


and  it  is  easy  to  establish  from  (B.16)  that 


m.  t    m+l,t-l   ^     '   ra.,  ^-l     m+l,  i-2 


(B.18) 


By  this  recurrence  relation^  wo  may  prove  from  (B.17)  by  induction  that. 


m,l      (m+P.+l);    ' 


(B.19) 


i.e. 


/  I   ^  •  '  '    m.'(m-hf.+l)  ; 


(B.20) 


r=l 
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Retiirning  to  (B.lj)  and  (B.I5)  we  find 


l?'^(-k,t)t2(k,s) 


/   I   m.'n.'(m+l+l).'(n+l+l)i 


U,-,(r  ,x   ) 
11^  o  o 


t  l+2m  l+2n 

(kv^t)    (kv^s) 


(B.21) 


2]  (^+1)^ 


11^  o'  o' 


'  -^i+i^^kv^^)  '^^+l(2l^l^) 


kv  t 


kv  s 


which  is  equivalent  to  (3.3l) 
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